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ABSTRACT 


A  semi-empirical  method  for  modelling  the  loss  of  electron  flaxes  in  the  earth's  mag¬ 
netosphere  was  developed.  An  equation  for  the  integral-energy  omnidirectional  electron  flux 
as  a  function  of  time  and  magnetic  Held  strength  was  derived  from  pitch-angle  diffusion 
theory. 


"’This  flux  equation  was  the  basis  for  a  computer  data-fltting  program  written  at  the 
Air  Force  Weapons  Laboratory  (AFWL)  to  fit  the  AFWL  Trapped  Electron  Da.a  flase.  The 
program  utilised  a  least-squares  fit  and  incorporated  random  variations  of  the  characteristic 
exponential  loss  times  about  their  initial  values.  Au  improved  table  of  initial  less  times  was 
compiled  for  use  with  the  program. 


'^Thc  derived  flux  model  showed  substantial  agreement  with  the  empirical  data  base. 
Representative  plots  of  computed  flux  over  raw  data  are  shown  for  lr values  of  2. it  ami  2.4. 


I.  INTRODUCTION 


Objective 

Modern  technology,  both  civilian  and  military,  depends  greatly  on  the  ose  of  satellite 
aystems.  It  is  thus  important  to  the  Department  of  Defense  to  understand  and  predict  the 
environment  in  which  satellites  may  operate,  in  order  to  determine  survivability  and  vul¬ 
nerability  requirements.  One  part  of  that  environment  is  the  electron  fluxes  which  may  exist 
in  the  trapping  regions  of  the  earth's  magnetosphere.  Such  fluxes  may  arise  from  natural 
sources  or  from  injection  by  high-altitude  nuclear  bursts.  Because  such  fluxes  have  great 
capability  for  inflicting  damage  to  sensitive  satellite  components,  even  at  great  distances 
from  a  burst,  understanding  their  behavior  is  vitally  important. 

There  are  a  large  number  of  theoretical  treatments  in  the  literature  which  describe 
various  electron  loss  mechanisms  from  the  earth's  magnetosphere.  However,  the  empirical 
calculations  of  actual  electron  losses  are  often  conflicting,  and  frequently  are  restricted  in 
their  areas  of  coverage. 

Since  the  cessation  of  high  altitude  nuclear  testing  in  1962,  loss  studies  have  been  limited 
to  electrons  injected  into  the  magnetosphere  by  natural  magnetic  disturbances  such  as  solar 
storms.  Thus,  much  of  the  theory  which  has  been  developed  since  1962  has  had  restricted 
opportunity  for  empirical  testing. 

The  objective  of  this  study  is  to  produce  an  improved  algorithm  to  calculate  losses 
of  electrons  from  the  magnetosphere.  The  algorithm  is  semi-empirical,  because  it  utilizes 
the  Air  Force  Weapons  Laboratory  (AFWL)  Trapped  Electron  Data  Base  (Ref.  37)  measure¬ 
ment  of  electron  fluxes  (following  the  old  high  altitude  nuclear  bursts)  to  determine  the 
characteristic  loss  times  of  those  electrons. 

The  loss  times  thus  determined  and  the  algorithm  for  flux  "decay”  will  then  be  incor¬ 
porated  by  AFWL  into  the  SPECTER  computer  code  (Ref.  7-9)  for  improved  flux  calcula¬ 
tions  of  satellite  environments. 


Scope 

This  study  is  limited  to  investigation  of  one  loss  mechanism,  pitch-angle  diffusion,  using 
a  theoretical  formalism  developed  principally  by  Shuts  (Ref.  47,46).  The  study  has  also  been 
necessarily  limited  to  those  regions  of  space  and  those  energies  covered  by  the  satellites  from 
which  the  data  base  was  complied.  The  loss  equations  thus  developed  in  this  study  are  only 
implicitly  functions  of  energy,  E,  and  L-coordinnte  (explained  in  Chapter  II). 

The  sheer  volume  of  data  and  the  limited  time  for  this  study  also  limited  the  testing 
of  the  algorithm  developed  to  a  representative  subset  of  the  entire  data  base. 


A  specific  solution  of  the  pitch- angle  diffusion  equation  is  derived  using  the  method 
proposed  by  Shuts  (Ref.  47).  This  solution  is  then  fitted  to  the  electron  flux  data  using 
a  least-squares  fitting  program  written  by  AFWL,  with  the  coordination  of  this  author. 
The  solutions  are  performed  for  integral  fluxes  (above  threshold  energy  of  the  satellite 
detector)  and  for  constant  L- value.  The  solutions  are  thus  explicitly  functions  of  magnetic 
field  strength,  B,  but  implicitly  functions  of  and  L  . 

Use  of  the  AFWL/NTCTS  computer  facilities  was  essential  to  this  study.  Well  over  72 
hours  of  actual  computer-processing  time  were  used,  and  over  1000  computer- produced  flux 
plots  were  examined. 

Presentation 

A  brief  Background  section  is  presented,  which  covers  some  physical  relations  and 
terminology  used  in  the  study  of  space  physics.  This  section  is  used  to  establish  uniform 
symbol  notation,  and  to  provide  a  common  starting  point  for  the  reader  unfamiliar  with  the 
specialised  language  of  space  physics.  It  may  be  skipped  by  the  reader  more  familiar  with 
the  topic. 

Chapter  m  outlines  the  theory  used  to  develop  the  flux  equations  which  were  fitted  to 

the  data. 

Following  the  theory,  chapter  is  a  chapter  outlining  the  method  of  data  analysis. 
Descriptions  are  given  of  the  data  base  and  of  the  AFWL  computer  codes  used  in  this  study. 

Representative  results  of  the  study  are  presented  in  Chapter  V,  in  both  tabular  and 
graphic  form.  Again,  the  sheer  volume  of  data  makes  inclusion  of  all  plots  and  fitted  values 
prohibitive. 

Conclusions  and  recommendations  for  improvements  for  further  study  are  presented  in 
Chapter  VI. 

The  Bibliography  includes  05  sources  and  provides  a  comprehensive  summary  of  existing 
literature  of  relevant  topics  about  radiation  trapped  in  the  magnetosphere  and  of  pitch-angle 
diffusion. 


E.  BACKGROUND 


The  equation  of  motion  of  a  charged  particle  in  magnetic,  electric,  and  external  gravita¬ 
tional  fields  is: 

d2? 


where 


m  dF  ~  mg  +  q 


df 


r—f(t,7 o,t»0);  —  =  5 

m  =  particle  mass, 

3  =5  gravitational  constant, 
q  =  charge, 

B  =  magnetic  field, 

£  =  total  electric  field 


(Ref.  11:23-24.) 


In  a  uniform  magnetic  field,  a  charged  particle  will  more  in  helical  fashion  along  and 
around  the  field  lines  with  a  cyclotron  radius,  pc,  (also  called  gyroradius  or  Larmor  radius) 
about  its  “guiding  center”.  The  particle’s  motion  may  thus  be  separated  into  its  rotation 
around  the  field  line  and  the  motion  of  its  center  of  rotation,  or  guiding  center,  along  the 
field  line  (Ref.  11:24-25). 

The  period  of  cyclotron  rotation,  rc,  is  defined  by  (Ref.  39:5): 

2jtto  2  repc 


T.  as 


qB 


v± 


where  m  is  the  relativistic  re-"**: 


fn  =  ffio7  =  mo 


and  where  »j_  is  component  of  velocity  perpendicular  to  B.  The  cyclotron  frequency  is 
simply  2v  over  the  period  (Ref.  39:5): 


_  2*  _  <75 
*  r,  m 


The  angle  between  a  particle’s  local  velocity,  5,  and  the  magnetic  field,  B,  is  its  pitch- 
ingle,  a,  which  is  defined  as: 


PA  •  ri) 

c*  =  arecosl  —  1 «  arcsin^ -----  j. 


>nn 


n  a  uniform  static  magnetic  field,  »;j  is  constant  and  i>j_  is  constant,  and  hence  n  is  constant 
rhe  particle  will  then  move  with  a  uniform  circular  motion  around  the  field  line  and  iinif< 
cctilincar  motion  along  the  field  line,  resulting  in  a  helical  motion  about  the  line,  since 

dp/dt  =  qd  X  If  =  m,(i, 

ind 

m3|)  =  q{v  X  5).|  =  0. 


0, 

constant, 

X  S)x, 
(qv±B)/m 


constant,  which  is  a  constant  centripetal  acceleration 
(Ref.39:4-7;  Ref.  11:23-38). 


The  concept  of  the  first  adiabatic  invariant  arises  naturally  from  the  gt  iding  center 
approximation.  The  first  adiabatic  invariant,  At,  (also  called  the  relativist  ic  magnetic 
moment)  is  defined  bj: 

A 

M  =  - - -  =  constant 

2moB 

where  jlj_  is  the  particle’s  perpendicular  momentum  in  the  guiding  center  ap  iroximation. 
The  assumptions  implicit  in  calling  M  invariant  are  that  the  spatial  variation  >f  li  is  small 
compared  to  pe  and  the  time  variation  of  B  is  small  compared  to  rc  (Ref.  39:1 1-23). 

If  a  particle  has  a  constant  velocity  along  a  field  line  in  the  guiding  center  frame  or 
reference,  then 

sin'a(s)  si  a2  a,- 

— ~  ~n —  =  instant 
B(s)  B{ 

where  s  is  the  arc  length  along  the  field  line  and  *  is  any  point  on  the  field  line.  This  assumes 
that  the  particle’s  kinetic  energy  remains  constant  as  it  follows  the  Geld  line,  if  a  particle 
enters  the  field  at  point  /?,-  with  pitch  angle  a;,  its  velocity  along  the  field  line  is 

«D (•»)  —  o  cos  a(e)  =  v  •  \J  1  -  “  -  sin2®,-. 

If  the  field  is  increasing  in  the  direction  of  V|,  then  the  “mirror  point"  of  the  particle,  whpre 
its  parallel  velocity  is  zero,  is 

n  ---  Bi  =  B{s) 

m  sin*a<  sin2a(s) 

and  the  particle  has  a  local  pitch  angle  of  909  at  that  mirror  point  (Ref.  39:34  -42). 

When  the  magnetic  field  has  a  geometry  like  that  of  a  dipole  field,  increasing  at  the 
poles  and  decreasing  in  the  midpoint,  then  the  particle  is  in  the  so-called  “magnetic  bottle”, 
and  is  trapped  between  the  mirror  points.  There  is  some  minimum  B  value  between  the 
mirror  points  which  is  called  Bo,  the  equatorial  field  strength.  Using  the  gliding  center 
approximation  and  approximating  the  earth’s  magnetic  field  by  a  dipole  fiel  I,  a  particle 
trapped  in  the  earth's  field  has  three  distinct  motions.  It  circles  rapidly  about field  line,  it 
bounces  along  the  line  between  mirror  points,  and  it  drifts  slowly  in  longitude  around  the 
earth.  All  three  motions  take  place  with  different  speeds,  so  they  arc  distinct  (Ref.  11:25; 
Ref.  29:34-85). 

From  the  mirror  point  definition,  and  from  the  velocity  equations,  we  sec  that 


and  that 


, .  .  fin*)  [m> 

u±(')  —  v  s.na,y  - --  =  ®  ^ 
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Note  that  the  particle’s  mirror  point  depends  only  on  its  "injection*  point  and  its 
injection  pitch  angle,  not  on  its  initial  energy  or  velocity.  A  mirror  point  is  a  consequence 
of  the  field  alone,  and  all  particles  injected  at  the  same  point  with  the  same  pitch  angle  will 
mirror  at  the  same  point  on  the  field  line.  Of  course,  this  is  only  true  if  no  external  forces 
are  acting. 


The  parallel  velocity,  »j,  of  the  particle  will  be  the  maximun  at  fi0.  wher?  /?{*)  is  the 
minimum.  Thus,  a  trapped  particle  spends  most  of  its  time  near  the  mirror  po  nts,  and  the 
least  amount  of  its  time  transiting  the  equatorial  Geld  regions  (Ref.  29:34-44). 

The  bounce  period,  »&,  of  a  trapped  particle  is  generally  much  greater  than  its  cyclotron 
period,  rt,  and  is  defined  by: 

„  ds  2  /*-»  d» 

n  =  2  —tt  —  - 1  — — — ~ 

J.mi  «,„(*)  vJ.mi  v/j TT(B(a)/Bm) 

where  and  «m,  are  the  mirror  points  on  the  field  line. 

As  a  particle  in  a  dipole  field  like  the  earth’s  bounces  along  the  field  line,  it  also  drifts 
perpendicularly  to  the  field  line,  due  to  external  forces,  field  gradients,  and  fie'd  curvature, 
ns  well  as  other  effects  such  as  time-dependent  field  changes.  This  drift  is  slew  compared 
to  the  bounce  period.  As  the  particle  bounces  and  drifts,  it  traces  out  a  surface  between  its 
mirror  points  and  around  the  earth,  called  a  drift  shell  (Ref.  39:9-19). 


The  concepts  of  the  second  and  third  adiabatic  invariants  arise  from  the  above  behavior. 
If  the  forces  acting  on  the  particle  remain  almost  constant  over  its  bounce  period,  the  iceond 
adiabatic  invariant,  J ,  is  defined  by 


Pi  * 


where  p|  is  the  momentum  component  parallel  to  B  and  d»  measures  arc  lenitth  from  the 
equator  (Ref.  45:12;  Ref.  (1:3-31).  If  the  forces  acting  on  a  particle  remain  almost  constant 
over  its  drift  period,  the  third  adiabatic  invariant,  <t>,  which  is  the  magnetic  flux  enclosed  by 
a  drift  shell,  is  defined  by 

«  =  j  A0  a  di 

where  *>  the  magnetic  vector  potential  and  where  the  integration  is  performed  over  a 
curve,  «,  which  lies  in  the  drift  shell.  The  third  invariant,  ♦,  is  defined  and  computed  for 
the  drift  shell  of  the  guiding  center,  with  constant  field,  and  not  for  a  drift  shell  which  a 
particle  may  physically  trace  out  under  short-term  conditions  (Ref.  45:12;  Ref  39:76-  79). 


The  guiding  drift  shell  of  a  particle  may  also  be  referred  to  as  an  "invariant  surface,” 
which  is  composed  of  field  lines  which  end  at  the  mirror  points.  The  three  adiaba.  ic  invariants 
uniquely  define  an  invariant  surface  (Ref.  6:3-31  to  3-32). 

The  more  common  set  of  parameters  used  to  define  a  particle’s  position  in  the  mag¬ 
netosphere  (or  to  define  an  invariant  surface)  is  the  B-b  coordinate  system.  The  fj  parameter 
is  magnetic  field  strength,  and  L  is  the  Mcllwain  b  parameter  defined  by  L  =  ?o //<k  where 
lip,  is  the  earth’s  radius  (s=s  8371  km)  (Ref.  8:3-33)  and  (Ref.  39:53). 


The  value  of  L,  therefore,  is  equal  to  the  distance,  in  earth  radii,  of  tl  e  equatorial 
point  on  a  field  line  in  a  dipole  field.  If  the  field  is  not  symmetric  (not  a  pc:  feet  dipole), 
the  invariant  surface  is  not  so  well  defined.  However,  for  most  field  lines  in  he  trapping 
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region  of  the  earth's  magnetosphere,  L  varies  by  less  than  one  percent  along  the  line,  so 
the  D-L  system  is  adequate.  The  lcss-than-exact  symmetry  of  the  earth’s  dipole  field  and 
outside  forces  may  cause  particle  drift  to  deviate  from  a  perfect  asimuthal  course.  This 
means  that  its  L  value  will  vary  in  the  course  of  a  drift  period.  This  variation  of  L  values 
is  not  significant  below  L  ==  3,  but  may  be  so  above  that  value.  The  average  //-value  of  all 
intersecting  invariant  surfaces  is  called  an  L-tkell.  This  Lr value  defines  a  set  of  surfaces  along 
field  lines,  which  may  end  at  differing  mirror  points.  The  //-shells  arc  considered  to  intersect 
the  earth’s  surface,  even  though  mirror  points  do  not  extend  through  the  at  mo:  phere  to  the 
surface  (Ref.  8:3-33  to  3-35;  Ref.  39:53). 

The  various  mechanisms  which  operate  to  cause  trapped  particles  to  change  //-shell  or 
to  be  lost  from  the  trapping  region  will  necessarily  violate  one  or  more  of  the  adiabatic 
invariants.  For  example,  pitch-angle  diffusion  violates  M  or  J ,  or  both;  and  radial  diffusion 
violates  (Ref.  45:48).  These  concepts  are  discussed  further  in  Chapter  III. 


m.  THEORY 


nr 


While  the  motion  of  particles  trapped  in  a  magnetic  field  is  welt  understood,  the 
mechanisms  of  losses  from  magnetic  fields  such  as  the  earth’s  magnetosphere  are  less  well 
understood.  It  is  generally  agreed,  however,  that  pitch-angle  diffusion  into  th<  loss  cone  is 
one  of  the  predominant  mechanisms  for  removal  of  charged  particles  for  mid-range  //-values. 
The  loss  cone  angle  is  the  lower  limit  of  pitch  angle  for  trapped  particles.  Any  particles  with 
smaller  pitch  angles  will  mirror  in  the  sensible  atmosphere  and  will  be  lost  by  atmospheric 
scattering. 

Roberts  (Ref.  41)  has  noted  that  a  pitch-angle  scattering  mechanism  must  be  extant 
for  pitch-angle  diffusion  to  occur.  Such  a  mechanism  would  necessarily  violate  one  or  more 
of  the  adiabatic  invariants. 


It  is  not  the  purpose  of  this  study  to  ascertain  the  true  source  mechani  im  of  pitch- 
angle  diffusion.  Such  a  study  is  beyond  the  scope  of  this  paper.  Indeed,  the  search  for 
this  mechanism  has  been  going  on  for  decades.  Roberts  (Ref.  41)  has  postulated  that  the 
“perturbation  forces”  causing  such  diffusion  may  result  from  turbulent  ambient  electric  or 
magnetic  fields  or  from  collisions  with  other  trapped  or  non-trapped  particles.  Lyons  (Refs. 
20-28)  has  formulated  extensive  theory  based  on  resonant  interactions  of  so-calh  d  “whistler” 
VLF  waves  with  the  trapped- particles’  gyrofrequencies:  “cyclotron-resonance* .  Regardless 
of  the  actual  physical  mechanism  behind  pitch-angle  diffusion,  its  general  treatment  is 
mathematically  the  same,  and  the  physical  results  are  the  same.  (Ref.  41:308)  This  study 
follows  the  methods  of  Roberts  (Ref.  41),  Shuls  and  Lanserotti  (Ref.  45)  and  Shuls  (Ref. 
47)  in  Misqunf  that  pitch-angle  diffusior  is  an  operative  process.  One  of  the  purposes  of 
this  study  is  the  formulation  of  the  equations  necessary  to  validate  that  assumption  against 
experimental  data,  and  to  perform  that  experimental  validation. 

It  should  be  emphasised  that,  while  pitch-angle  diffusion  is  assumed  in  this  paper  to 
be  the  dominant  loss  mechanism  for  b- values  and  altitudes  considered,  it  is  not  the  only 
loss  mechanism.  Several  investigators  have  formulated  radial  diffusion  (cross^-L)  theories 
(for  example,  Walt  (Ref.  58)  or  Tomassian  (Ref.  54)),  coulombic  or  collision-scattering- 
diffusion  theories  (for  example,  Wentworth  (Ref.  80))  and  multiple  diffusion  theories,  such  as 
combinations  of  pitch-angle,  energy,  and/or  radial  diffusion  (for  example,  Walt  iRef.  58)  and 
Lyons  (Ref.  20)),  all  of  which  show  some  agreement  with  experimental  data,  la  particular, 
atmospheric  scattering  is  obviously  a  dominant  force  at  very  low  L  values  (Ref.  30). 

The  introduction  of  multiple  low-mechanisms  makes  explicit  solutions  of  iny  diffusion 
equation  extremely  difficult.  To  simplify  the  problem  and  render  it  amenable  to  the  method 
of  Shuls  (Ref.  47),  this  study  assumes  that  “as  a  rule,  radial  diffusion  enables  the  radiation 
belts  to  become  populated  from  an  external  source  (or  rearranges  particles  injected  by  an 
internal  source),  while  pitch  angle  diffusion  causes  particle  loss  to  an  atmospheric  sink”  (Ref. 
45:48). 

Pitch-angle  diffusion,  while  simple  in  concept,  is  complicated  in  detail.  Roberts’  (Ref. 
41:307-337)  treatment  of  the  general  mechanism  is  particularly  descriptive: 


At  the  magnetic  equator,  a  particle’s  pitch  angle,  ao,  is  determined  by 


2  =  COS  O.0 


_ Pi.  _ 
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-where  and  Pj_f  are  the  parallel  and  perpendieluar  components  of  the  particle's 
momentum  with  respect  to  the  field  line  at  the  equatorial  value  of  fio-  Roberts  next  defines 
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where  mo  is  the  rest  mass.  These  are  defined  since  the  first  adiabatic  invariant,  M,  is 
proportional  to 

(3)  M~\ 

and,  if  Wj,  <  W±,,  the  second  adiabatic  invariant,  J,  is  approximately  prrportional  to 

WW. 

rDWx, 
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where 


ra  =  the  particle’s  bounce  period, 


relativistic  mass  m 


rest  mass 


A  particle’s  path  may  be  defined  by  plotting  its  values  of  Wj_t  and  W|,  as  in  Figure  1. 


W||0*  P*no  /2mo 

Fig.  1.  Loss  of  s  trapped  particle  by  a  random  walk  into  the  pitch- 

anale  loss  cone.  (Ref.  41:306) 


Note  that  lines  of  constant  equatorial  pitch  angle  pass  through  the  origin,  because 


Any  force  which  violates  either  M  or  /,  the  first  or  second  adiabatic  invariant,  will 
produce  ‘diffusion*  of  the  particle’s  pitch  angle  as  shown  in  Figure  1  as  Wj_,  and  W|( 
change  with  time.  When  a  particle's  equatorial  Ditch  angle  diffuses  to  the  valie  of  the  loss 
cone  angle, 

(^)  * t  ”  cos  a,  **  \ll  - 

V  c*ioo 

where  Btoo  is  the  value  at  100  km  altitude,  then  it  is  lost  by  atmospheric  scattering.  "The 
loss  cone  serves  as  a  ‘sink’  for  particles*  undergoing  pitch  angle  diffusion  (Ref.  41:307). 

Roberts  makes  two  very  important  points  about  pitch-angle  diffusion  as  a  oss-mechan- 

ism: 


(l)  with  no  source,  the  entire  radiation-belt  would  be  depleted  of  particlcr, 

(3)  the  loss  cone  approximation  is  just  that,  since  the  atmosphere  is  not  sharply  defined 
at  100  km.  However,  if  atmotpherie  pitch-angle  scattering  is  not  the  primary  focus, 
the  loss  come  approximation  may  be  useful.  ‘Naturally,  when  the  loss-:one  concept 
is  used,  detailed  agreement  between  theory  and  experiment  cannot  be  expected  in 
the  region  near  the  edge  of  the  loss  cone*  (Ref.  41:308). 

Both  Roberts  (Ref.  41)  and  Shuls  and  Lanzerotti  (Ref.  45)  develop  the  pitch  angle 
diffusion  equation  from  the  Fokker-Planck  equation  using  a  particle  distribution  function, 
?• 


The  Fokker-Planck  equation  is  a  formalism  which  arose  from  the  study  of  Brownian 
motion,  and  which  is  used  frequently  in  transport  theory  (Ref.  43:308).  The  characteristic 
Fokker-Planck  equation  for  trapped  particles  (which  ignores  radial  diffusion)  h 


eJ  =  _±o_ 

dt  7 p  dE 


+ 


i  e_ 

X  T[y)  9x 


*  T(y)  D„ 


—  — 

7 p  9E 


'IpDsB 


d l 

dE 


where  the  first  term,  subscripted  v,  represents  non-stochastic  (mean)  energy  loss  to  the 
atmosphere,  the  second  term  represents  pitch-angle  diffusion,  and  the  third  ter  n  represents 
range-straggling  (energy  diffusion).  The  term  T(jr)  is  defined  below,  and  the  Due  and  Dzt 
terms  are  the  characteristic  diffusion  coefficients;  7  is  relativistic  mass  ratio  (Ri  f.  45:55-58). 

Roberts  uses  a  distribution  of  particles  in  a  “tube*  of  force  about  a  0.  Id  line,  and 
Shuts  and  Lanserotti  use  a  phase-space  density  attribution  function  which  s  essentially 
equivalent  to  Roberts’.  Since  this  study  follows  the  methods  of  Shuls  (Ref.  4")  and  Shuls 
and  Lanserotti,  a  discussion  of  phase  space  is  necessary  (Ref.  45:15-23). 


Any  moving  particle  can  be  described  by  specifying  its  three  position  coordinates  and 
its  three  canonical  momentum  components.  This  completes  a  six-dimension  “phase-space” 
in  which  a  particle  can  travel  in  time. 


If  there  exists  a  system  of  a  large  number  of  particles  in  phase-space,  th«  system  can 
be  described  by  a  six-dimensional  distribution  function  where 

Pi  (»  **  1,2,3)  are  canonical  momentum  components, 

7,  (t  =  1,2,3)  are  position  coordinates,  and  t  is  time. 
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Thru,  /  is  the  number  of  particles  instantaneously  occupying  5-D  volume 

dPPd*^.  According  to  Liouville’s  theorem,  the  phase-space  volume  containing  the  system 
of  particles  moves  incompressibly  through  phase-space. 

Since  P  is  an  awkward  quantity  to  deal  with  physically,  Shulz  and  Lanzerctti  note  that 
P  transforms  to  the  more  familiar  p  as 


where  A  is  the  electromagnetic  vector  potential. 

Hence,  f(P,<i,t)  —  /(?,?,<)  since  <J  =  ?  and  the  P  to  p  transformation  has  a  unit 
Jacobian.  The  position-momentum  distribution  function  /(?,  ?,  t)  defines  the  particles  oc¬ 
cupying  the  8-D  volume  d?pds?  at  any  given  t.  The  pitch  angle  diffusion  equation  is  given  in 
terms  of  a  phase-space  distribution  function  which  is  numerically  equal  to  the  more  easily- 
definable  position- momentum  distribution  function  (Ref.  45:15-22).  At  constant  energy  and 
£-shell,  and  under  the  action  of  some  source  S,  the  diffusion  equation  in  phase-  <pace  can  be 
written 

rc  ft=7mUxT[y)D^hs 

where  T(y)  is  the  quarter- bounce  integral  path  length  function  of  vT  -  z2  where 
i  a  cosoi-  The  quarter  bounce  integral  gives  the  length  of  the  trajectory  of  a  particle  (in 
units  of  JZb,  the  distance  to  the  equatorial  crossing)  from  the  equator  to  the  mirror  point. 
The  exact  definition  of  T  is 

r./  — 

Jo  cos  a 

where  s  is  the  distance  along  the  field  line  and  a  is  the  local  pitch  angle  (Ref.  10:4020-4030). 

D„  is  the  bounce-averaged  diffusion  coefficient  as  a  function  of  z  (Ref.  41  and  Ref. 
45).  The  boundary  conditions  are  that  J  —  0  at  z  =  x,  (some  cutoff  value)  and  J{z  =  0)  is 
finite. 

Shuls  and  Lanserotti  (Ref.  45:182)  rewrite  the  above  as 


Oj  i 

dt**  x  dz  l* 


Z 

-  ~D, 


and  further  approximate  that  the  second  term  is  negligible  for  z  <  1  since  z2  +  y2  1. 
If  T[y)  «  T(l),  the  second  term  disappears  and  the  equation  is  a  diffusion  equation  in 
cylindrical  coordinates.  (These  approximations  are  not  used  in  this  work.  This  development 
is  used  to  aid  in  understanding  the  exact  solutions  which  follow  from  Shuls  Ref.  47)).  If 
Da,  and  S  are  then  independent  of  s,  and  if  J(x,t)  =  X(x)  T(t),  the  eigenfurctions  of  (8) 
arc  Bessel  functions  of  order  zero.  The  general  solution  to  (8)  would  then  be 


?(*»*)  =*  /»(*)  +  £«»(*)  Jo 


where  7<»  >>  the  steady-state  solution 


Mast -(0} 


10 


The  Kn  are  the  zeroes  of  Jo  (n  =vl,  2, . . .)  and  the  o»(t)  vary  as  e~( •/*•),  where  r«  =« 
z3/(Z)3<e3).  Thus,  the  pitch  angle  distribution  (and  hence  the  directional  or  omnidirectional 
flux)  is  shown  to  be  the  sum  of  a  steady  state  and  higher  order  eigenmodes  (Ref.  45:160-168). 

The  steady  state  can  be  thought  of  as  the  "normal”  or  “quiet-time*  value  of  flux  which 
exists  in  equilibrium  with  the  source  S.  An  "injection”  of  particles  would  then  result  in  a 
perturbed  distribution  function  with  several  eigenmodes,  each  with  a  characteristic  decay¬ 
time  r„.  The  higher  order  modes  decay  faster,  and  eventually  only  the  fundamental  mode 
would  remain,  which  would  decay  exponentially  to  reach  the  steady  state  (after  "infinite* 
time)  (Ref.  41  and  45).  There  is  ample  evidence  that  such  a  process  does  indeed  occur. 
Roberts  (Ref.  41)  cites  Explorer  XV  data  for  the  28  October  1082  Russian  explosion  as  an 
example.  Rosen  and  Sanders  (Ref.  43)  also  note  that  decay  is  faster  immediately  after  solar 
magnetic  storm  activity  than  during  quieter  periods.  The  similarities  in  Figures  2  and  3 
illustrate  the  way  in  which  a  decay  of  higher  eigenmodes  (Figure  2)  can  approximate  the 
temporal  evaluation  of  an  actual  electron  distribution  (Figure  3). 

The  primary  difficulties  with  the  above  discussion  arise  from  the  approximations  that: 

(1)  T(y)  =  T(l)  =  constant 

(2)  z3  <  1  (orytm  1) 

(3)  Dxx  varies  little  (or  not  at  all)  with  z. 

Roberts  (Ref.  41:311)  uses  the  argument  that  the  full  spiral  path  (and  hence  T{y})  varies 
only  by  a  factor  of  1.4  as  z  varies  from  0  to  0.9  and  by  1.9  as  z  varies  from  0  to  1.  Shuls 
and  Lanserotti  (Ref.  45:163)  use  the  simplifying  assumptions  that  T(y)  ==  T(l)  and  Dxu  is 
independent  of  z.  Both  use  the  assumption  that  z3  <  1. 

The  present  study  requires  that  z  be  allowed  to  vary  from  0  to  zc  (up  1o  near  1)  in 
order  to  adequately  examine  the  AFWL  trapped  electron  data  base  (Ref.  37).  In  addition, 
if  z  is  not  small,  the  approximation  that  T[y)  =  T(l)  is  a  poor  one,  since  y  is  not  close  to 
1. 


The  function  T[y)  can  be  shown,  within  0.57 %  to  be  approximated  as  (Rc  f.  10:4030) 
(ll)  T(y)  «  T(0)  -  (r(0)  -  r(l)]y5/4 

where  7*(0)  =  1 J-  |(ln(2  +  V5))/(2>/5)j  «  1.3802  and  T(l)  =  (jrv^/B  «  0.7405  (Ref. 
45:19). 

Obviously,  if  z  approaches  1,  y  approaches  0  and  the  T{y)  =  T(l)  approximation 
becomes  invalid. 

The  three  limiting  approximations  are  removed  in  the  treatment  of  Shulz  (Ref.  47)  by 
introduction  of  a  new  "canonical”  variable  z  such  that: 


(12) 


Z(v)  =  /  *'  7V)  dz'  *  /  /r(yW 

JO  Jf 


and  a  corresponding  diffusion  coefficient 


Figure  2.  Postulated  Decay  cf  Higher  Flux  Eigenmodes  to  the  Steady  Strtc  by  Pitch- 
Angle  Diffusion  (Ref.  41:313). 
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Figure  3.  Temporal  evolution  of  1.9  Mot  Electrons  from  the  Russian  3  Hurst  (Ref- 


Using  (li)  in  (12)  sad  performing  the  integration: 


(14) 

where  7( 0) 


z  « 


^l[r(oy-(r(o)-T(i)]//4 
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10/35  and  7(1)  —  0  (Ref.  46:5213). 


Using  (7),  (12)  and  (13),  the  diffusion  equation  becomes: 


(15) 


s 2 

at 


oj 

as 


+  s. 


Shnl*  (Ref.  47:0)  states  that  (15)  is  a  “canonical”  diffusion  equation  in  that  there  is  no 
Jacobian  factor  which  fails  to  commote  with  DMM  as  was  the  case  in  (7). 

The  assumptions  which  now  must  be  made  are: 

(1)  Z>„  is  some  “suitably  simple  function  of  z”  (Ref.  47:0)  so  that  exact  ei  jenf unctions 
f„(z)  may  be  specified, 

(2)  S,  the  distributed  natural  environment  source,  is  independent  of  z. 

Assumption  (1)  is  the  limiting  assumption,  since  if  the  function  of  z  is  very  complex,  the 
diffusion  equation  becomes  extremely  difficult  to  solve.  Assumption  (2)  is  both  simple  and 
reasonable,  however.  The  source  of  the  natural  environment  must  be  close  tc  constant  in 
order  to  be  the  “driver”  of  quiet-time  equilibrium.  It  must  also  be  distributed  'airly  closely 
to  the  steady-state  quiet  time  distribution  for  the  same  reason. 

Shuls  (Ref.  47:7)  states  that  even  if  D,a  is  not  of  a  functional  form  to  yield  exact 
eigenfunctions,  there  may  be  a  t>„  which  resembles  D„  "closely  enough”  and  for  which 
exact  eigenfunctions  f„(z)  (resembling  y„(z))  are  known. 

The  following  derivatior  of  the  exact  eigenfunctions  and  the  corresponds  x  omnidirec¬ 
tional  differential-energy  flux  arises  directly  from  Shuls  (Ref.  47). 


The  first  assumption  by  Shuls  is  that 


where  o  is  some  number,  not  necessarily  an  integer,  less  than  2,  and  where  Dt,  is  the  value 
of  D„  at  some  z«  <  10/35  where  /  vanishes.  Rol>erts  (1969)  also  uses  this  as  wrapt  ion. 

This  form  of  Dtt  allows  a  basis  set  of  orthogonal  eigenfunctions  to  be  si  own  for  the 
interval  of  interest:  0  <  z  <  ze.  Following  the  notation  above,  we  seek  st  me  D„  and 


I  %* 


eigenfunctions  gn[z)  which  "closely  resemble"  the  true  functions  and  which  satisfy 

<17) 


with  the  substitution 


5n(*)  =  Zaw{Pz() 


then  (17)  becomes  Bessel’s  equation.  The  exact  eigenfunctions  of  (17)  are  given  by 


'■"-UM  4-(i)  ) 


where  Km  are  the  nth  seroes  of  the  Bessel  function  /„  of  the  order  v,  where 


The  given  by  (18)  are  normalised  so  that 


f  9n  9m  dz  —  6„ 

J  o 


The  eigenvalues  of  (18)  arc  given  by 


«r)2D„.. 


Shuts  and  Boucher  (Ref.  18:6)  state  that,  since  particles  are  not  lost  at  r  =  0  (but 
rather  at  t  =  z„),  the  diffusion  current  must  go  to  sero  at  *  =  0.  A  diffusion  current  at 
g  mt  0  would  imply  diffusion  of  pitch  angles  (and  hence  mirror  points)  into  the  equator, 
which  is  not  reasonable.  The  condition  of  no  diffusion  current  at  z  =  0  corresponds  to 
lims_4)  D',,$'n(z)  —  0. 

Shuts  and  Bonchcr  (Ref.  18)  ose  the  series  expansion  of  /*  from  Abraraowic ;  and  Stcgun 
(Ref.  I)  to  show  the  limit  as  z  approaches  sero.  This  author  shows  the  limit  directly  for  the 
special  case  eigenfunction  which  is  used  in  this  report.  The  derivatior  is  shown  in 

Appendix  A. 
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In  order  to  arrive  at  the  differential-energy  omnidirectional  flax,  Shall  (Ref.  47)  assumes 
that  /(*,  t)-»/oo(*)  as  t— »oo;  that  is,  the  solution  distribution  function  has  some  steady-state 
value.  This  also  assumes  that  $  is  constant  in  time,  but  not  in  x  or  z.  The  distribution 
is  necessary  in  order  to  expand  7  as  a  scries  of  eigenfunctions  in  the  flux  equation.  To  arrive 
at  an  expression  for  7oe  start  by  setting  dj/dt  to  zero  in  (15)  and  integrating  twice, 
which  gives: 

<20) 

Expanding  S  as  a  series  of  orthogonal  eigenfunctions  gives: 

(20a)  $(*")=  f>„5„(-~") 

»*0 

and  using  orthogonality  (18a)  with  (20a)  gives: 

(20b)  a.  =  /  'S(^)  gjz”)  dz>". 

J  o 

Substituting  (20a)  into  (20)  gives 

(  2°C  )  /„(*)  .  I'-jfcl'  £  U*)  **  <*■ 

Rearranging  the  order  of  summation  and  integration  in  (20c)  gives: 

(20d)  ?<»(*)  =£>»/  D^fo 

Now  integrating  (17a)  twice  gives: 

M*l-  C^L'^*** 


or  since  is  a  constant, 

<2°e)  v-f&jcw-- 

Substitution  of  (20e)  into  (20d)  gives: 

(20f)  7.w- 

n-0  A" 


Substitution  of  (20b)  into  (20f)  gives: 

(  21  )  ?»(*)  =  £  W  fV)  Intf)  dzf 

n-0  X" 

where  S(^)  is  positive  (hence  a  source). 

Shul*  and  Lanzerotti  (Ref.  45:39-40)  point  out  that  the  directional  flux  = 

Pi  /O’IjPX*^)  whero  /  is  the  same  distribution  function  as  in  this  paper.  8h  ill  (Ref.  47) 


points  out  that  j  most  be  integrated  over  all  angles  in  the  unit  sphere  in  p-spr  ce  to  obtain 
the  omnidirectional  flux  J: 


/•cot  a, 

J  =  4  wp2  I  7  d(cos  <*)• 
Jo 


Since 

sin2a  sinsao 

B  ^  Bo 

and 

sin2a(§), 

it  follows  that 

and  that 

s/l-y*(B/Do) 

Now  since 

if  eosa  =*  cos  a,, 

y*  =  y«2;  if  cos  a  *  0,  y2  =  Bo/B  and  (22)  becomes 

=i|  .r*  -mwfdtf) 

P  JIB./B )  J\  —  U*(B/Iin) 


which,  upon  combining  terms  and  reversing  the  limits  (•): 

(23)  ,  n  2  f{B*,B)W  Do)?  djv2) 

A.-  y/l-vHD/Bo) 


Vs(B/Bo) 


Integrating  (23)  by  parts 


judv  =  ut?  -  Jvdu  where  u  =*  J 
•  _  (B/Bp)  d(y2) 


v  =  -2  \/i-ir 


1-vHB/Bc)' 
L  J  B\ 


*C  V' “''(£)  w)*""5 


and  since 


df  Of  _dj_ 

d(y2)  Dz  ,%2) 


'■  ] 


a 


s  ** 


h-. 


>  ■ 


i 


'will  give  (since  the  definite  integral  evaluates  to  tero,  because  7  is  0  at  ye  by  definition): 

But,  from  (12),  dz  =  —y  T(y)  dy.  So 

_  -V  T(y)  dy  _  -yT{y)  dy  _  — T*(u) 

<%*)  d{y-)  2y  dy  2 

Hence,  (24)  ultimately  becomes: 

(25) 


Shulz  then  expands  7  as  a  weighted  series  of  eigenfunctions: 

OO 

(2«)  7  =  /«(*)  +  £  >U(£,  i; «)  ?„(*)• 


ft— 0 


This  allows  v  he  flux  7  in  (25)  to  be  written: 
(27) 


J  =»  2rp2 


/•<»•/»)  / 7n\ 


-  f;  a.(b,  L-,  u»>  nw 


The  following  two  expressions  are  derived  in  Appendix  B. 

(28)  «o  ^  /■*. 

/«(*)  =  _ 

n-0  An 

v-(*+D/2 


/«(*)  -  £  ®  /*$  5„(*')  d(/)t 

_ A  Af.  »0 


(29) 


fni*)  ~ 


y/{2  -  ff)/ze 


,^,:(KX5r"'-Nin 

.ar^nw-ifeni 


The  actual  form  of  the  flux  equation  awaits  the  choice  of  o,  which  also  determines 
the  dependence  of  Dtt  on  z.  Shulz  uses  a  =  0  for  purposes  of  comparison  of  several 
computational  methods.  This  choice  is  also  attractive  because  it  makes  D,t  independent 
of  z  (see  17(b))  which  was  a  major  assumption  to  begin  with: 


(30) 


Dtg  —  D**,  • 


This  choice  is  made  in  the  present  study  for  the  above  reasons  and  also  for  the  sake  of 
simplicity  in  computations. 
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With  tr  =  0,  the  Bessel  functios  within  eigenfunctions  become  order  —1/2,  which  can 
be  simply  written  as:  _ 

=  y[~x  cos(R), 


(31) 

(32) 


*R'-\/2/(xR) 


and  the  seroes  k„  occur  at  (2n  +  1)(t/2). 


This  simplifies  (32),  since  the  second  term  will  always  be  sero  if  /'_l/2(«n)  is  evaluated 
at  k„  =  (2n  +  t)(Jr/2). 


Equation  (32)  becomes 
J\ 


-1/2  (*") - ^ 


x(2  n  •+•  1)(jt/2) 


sin((2n  +  l)£) 


(33) 


'-*(*")“  \] *(2n  + 


mm 


(-i) 


n+l 


The  integral  with  respect  to  /  within  (28)  would  thus  be  evaluated  ar: 


Jo  -vfeWsi  sin((2n+i)?)  v  c  } 


•L 


*•  V^/ze  cos  ((2n  +  l)(>r/2se)  z’)dzl 
o  -sin  ((2n  +  l)(>r/2j) 

y/2/zc  sia((2n  +  l)(x/2)(z/ze)) 

((2n  +  l)(s-/2ae))(— l)"+l 

__  y/ljzc  sin((2n  +  l)(w/2)) 

”  ((2n+  l)(sr/2ae))(— l)n+r" 

The  5v„  in  (28)  would  be  evaluated  as  (from  (19)): 

...  , 

(35)  X„  =  (2n+1 


(34) 


Ki)  d,m" 


The  flux  equation  (27),  using  (31),  (33),  (34),  and  (35),  will  then  simplify  to  the  following 
for  ff  =  0: 


^ABt/S)  3  Vl -/*•)  «in((2n+l)(»/2)) 

_  /  V'°°  Ign+OU/g*, )(-»)’ fl 

4?  Z-"-°  (2«+  l)-(,r/2^)2  U„, 


\J irin+liTV/gi  ^ 


{(iXi)  ”•('*• - °i ;) 


+  £  (  \/*(2n  +  l)(jr/2){*/ee)  am((2n  ' '  2  ^ 


_ cos((2n+  t  )(x/2)(z/ze )) _ X 

*((2n+l)(*/2)Wze))2 

™  *■>] 

JfA  wi,. _j3!Z3_ 
LtS  *'  V27W2^T)|W)  (-D' 


2 Zc\zcJ  y  l)(»/2)(*/*«) 


_ V^) 

vv(*(*n+i){*M  (-ir+' 


— — —  cosf(2n  +  l)-*-*-) 

2)(*/*e)  V  %  zcJ 


\fr,  (-»"(<2” 


+  l)2zJV  r(2n  +  l){w/2){z/Sc) 


_ cog((2n  +  l)(y/2)(g/ 

)r((2n  +  l)(*/2)(z/ze)f  yjf/jxt 


•  l)(*/2)(*/z.))  A 

:e))2v/2/(T(2n  +  I  )(*/*))  J 


*  ™  <«  ]}• 

Equation  (38)  reduces  to  the  following  upon  simplification  and  cancellation  of  terms: 

"HirS'KSSSSW  •  \RE  »<* 
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The  first  term  in  Equation  (37)  represents  a  steady  state  value  of  the  flux,  ar  d  the  second 
term  represents  a  sum  of  eigenmodes.  Each  An  (Kef.  45:162)  has  the  form 
As  n  increases,  the  rn’s  become  shorter;  i.e.,  the  higher  eignmodes  decay  fatter,  so  that 
eventually  only  the  fundamental  mode  is  left,  which  decays  exponentially.  At  t  =  oo,  the 
fundamental  mode  has  reached  sero,  and  only  the  flux  remaining  is  the  steady  state  flux, 
which  corresponds  to  the  natural  environment  background. 

The  remaining  substitutions  in  (37)  are  ze,  z,  T(y),  and  ye: 

Vc  =  J - * - — , 

ViV« -(*/*>) 

(38)  =  =  Z(y)  =  ^  (1  -  y2)T(0)  -  ^(T(0)  -  T(t)](l  -  y<n/4>), 

ze  =  Z(Vc)  =  l-(l  -  y*2)T(0)  -  ~ [7*(0)  -  T(L)l(l  -  ye<"'4>) 

(as  shown  in  Equation  (14))  T(y)  is  given  by  Equation  (1 1). 

The  variable  of  integration  in  (37)  is  taken  as  y  rather  than  y2  for  -implicity  of 
integration. 

**  ie  preceeding  derivation  has  explicitly  followed  Shnls  (1981).  However,  Shuts  did 
not  slow  how  the  differential-energy  flux  related  to  the  integral-energy  flux.  Shuls  and 
Lanserotti  (Ref.  45:163)  stale  that  the  integral  flux  "will  then  scale  as*  the  differential  flux, 
without  derivation.  Since  the  AFWL  trapped  electron  data  bate  (Ref.  37)  exists  primarily  as 
omnidirectional  fluxes  integrated  over  broad  energy  bandpass,  the  relation  must  be  derived. 

It  will  be  shown  in  the  following  and  Lanserotti  (Ref.  45)  are  indeed  :orrect:  The 
differential-energy  flux  equation  (37)  will  scale  as  the  integral  flux,  and  the  energj  dependence 
will  be  imbedded  in  the  constants  Tor  the  source  and  the  eigenmodes. 

From  Roedercr  (Ref.  39:88),  the  integral  energy  omnidirectional  flux  can  l»c  written: 

J>E-  [  JdE 
Je 

where  /  is  the  differential-energy  flux,  as  in  the  previous  equations,  and  J>e  is  the  integral 
energy  flux  above  threshold  energy  E. 

A  principal  assumption  will  be  that  energy  and  pitch  angle  are  both  independent 
variables  of  the  flux.  Thus,  (26)  would  really  look  like 


(M  »)  ?(*,  E)  =  ?«,(£)  J^z)  +  £  An(E,  L\  t)  j?B(*) 

n«»0 

and  hence  in  (25)  dj/dz  would  look  like 

(39  b )  al  =  E )  7U(=)  +  £  ME,  h,  t )  U-). 

n-0 

Using  (27),  (39a)  and  (39b),  Equation  (39)  becomes,  at  a  given  L-shell: 
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'  — ; 

J>B  »  -2r£° pVoo,^) ^  Jmi ’  7U,(*)  /£°°  P2  £’  /„(*) * 

where  Toe,(x)  an<^  represent  all  of  the  angular-dependent  functions  of  z,  r,  or  z  shown 
in  (25)  through  (27). 

Using  the  previously  stated  assumption  (Ref.  45:152-103)  that  AnL{E)  arc  of  the  form 
an{E)t~,^T^h:^ ,  40)  becomes: 

J>  B  =  -2a  f  P'JootiE)  dE  f  ’  7oo,(*)  dx 

J  B  J*  i 

(41)  P2*n{E)*-,/,T'(E))  dE  f  '?n(z)dz 

a-o'* 

at  a  given  //-shell. 

To  "remove”  the  energy  dependence  of  the  decay-time  terms,  the  expectation  value 
of  e_|/lT»l  is  computed,  assuming  that  the  second  term  in  (41)  represents  a  listribution- 
f unction  of  the  energies. 

m  V  /  ST  P»«.(B)CK 


Then  (41)  can  be  written,  using  (42): 


J>B  «  2x[-  fK  p'J^E)  dE ]  jH  7UW  dx 
(«)  -2>rE  [/g  p2aB(£)  (e',/,r'’)t  / '  j'J*)  <**• 


Now  since  it  is  assumed  that  the  energy  and  angular  shape  functions  are  separable, 
the  / «*,,(£)  must  include  only  those  parts  of  (28)  which  have  energy  dependence.  The  only 
terms  in  (28)  which  could  have  any  energy  dependence  are  B  and  X„. 

Shuls  (Ref.  47:23)  notes  that,  while  7oo(*)  resembles  go(z)  in  functional  form,  7<»(~) 
will  coincide  exactly  with  ga(z)  only  if  S  is  directly  proportional  to  ju(z).  Fcr  simplicity, 
the  assumption  is  made  in  this  study  that  S  is  constant  over  the  interval  0  <  z  <  zc.  In 
physical  terms,  this  really  corresponds  to  a  steady-state  isotropic  source.  It  is  known  that 
the  natural  environment  is  not  isotropic  over  cither  space  or  lime,  but  the  variations  are 
generally  small  compared  to  the  variations  from  an  injection  by  a  nuclear  evert. 

Now  if  S  and  X„  are  the  only  terms  in  (28)  which  may  have  energy  dependence,  then 
7oe|(g)  in  (43)  would  correspond  to 

and  since  the  only  portion  of  the  which  may  have  any  energy  dependence  I'rom  (10))  is 
then  7^ ,(/•;)  should  be  proportional  to 

3(E) 
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where  Dt0,(E)  may  be  a  constant  in  E.  There  may,  of  coarse,  be  some  other  functional 
energy-dependence  included  in  /<*,,(/?)  which  is  yet- undetermined. 

Equation  (43)  is  generally  of  the  tame  functional  form  as  equation  (37)  if  the  terms  in 
square  brackets  in  (43)  correspond  to  the  combined  constant  terms  in  (37)  as  demonstrated 
above  and  as  shown  in  Appendix  D.  Equation  (43)  also  bears  out  Shuts  and  Lanzerolti's  (Ref. 
45:183)  prediction  that  the  integral  (lux  scales  with  the  differential  flux;  i.e.,  the  ‘pitch-angle 
shape*  dependence  is  invariant  for  differential-  or  integral-energy  fluxes. 

Since  it  is  impractical  to  evaluate  an  infinite  number  of  terms  in  (37),  only  the  terms 
up  to  n  —  3  will  be  shown.  Shulz  and  Lanzerotti  (Ref.  45:183)  state  that  the  I  igher  modes 
vanish  for  n  >  2.  This  in  fact  means  that  the  higher  modes  decay  so  rapidly  that  modes 
higher  that  »  =  3  should  not  be  seen  except  at  very  early  times.  The  final  form  of  Equation 
(43),  taken  to  n  =  3,  is  shown  in  Appendix  C.  This  is  the  form  of  the  flux  eqm.tion  used  in 
the  data  analysis,  except  that  fewer  terms  were  used  than  are  shown  in  Apperdix  C. 

The  general  solution  for  the  differential-energy  flux  (which,  from  (43),  is  directly  propor¬ 
tional  to  the  integral-energy  flux),  using  Bessel  functions  of  arbitrary  order  and  with  it- 
dependence  explicitly  shown,  is  given  in  Appendix  D. 


IV.  ANALYSIS  OF  DATA 


The  AFWL  Trapped  Electron  Date  Bate  consists  of  a  act  of  computer  tapes  of  aatellite 
counting  data  which  hare  been  assembled  from  hiatorieal  arc  hire*  (such  as  i  be  National 
Space  Sciences  Data  Center,  NSSDC)  (Ref.  37).  The  data  cover  injections  of  el  -ctrons  from 
nine  high- altitude  nuclear  detonations  shown  in  Table  I  (Ref.  0:8-2). 

The  Data  Base  has  been  organised  by  Pfitser  into  a  coherent  set  of  tapes  which 
present,  by  satellite,  net  electron  omnidirectional  integral-energy  (above  threshold)  fluxes  as 
a  function  of  time,  B,  and  L-sheil.  Tbc  computed  errors  for  each  data  point,  in  units  of  flux, 
are  also  given.  The  background  which  was  subtracted  to  give  net  flux  is  also  given  (Ref. 
37:11,  188-189). 

Additional  data  have  been  collected  by  AFWL  from  plots  published  in  early  literature. 
These  plots  have  been  photographically  enlarged,  and  the  data  points  digitised  and  placed 
in  computer  files  with  formats  similar  to  the  Pfitser  data  (Ref.  12:870-  871,  Ref.  31:848-483, 
Ref.  38:637-838). 

The  satellites  which  provided  the  data,  and  their  orbital  parameters,  are  sh  rwn  in  Table 
II  (Ref.  44:41-45). 

The  analysis  of  the  data  base  required  development  of  three  computer  programs  by 
personnel  at  AFWL,  in  coordination  with  this  author: 

(1)  "Program  DTABASE"  to  read  raw  data,  organise  it  for  processing,  and  output  it 
into  uniformly  formatted  files  (Ref.  15). 

(2)  A  program  to  plot  the  raw  data  and  to  plot  flux  curves  generated  by  the  fitting 
program. 

(3)  "Program  Elcctrofit",  a  fitting  program  to  take  data  points  from  DT  XBASE  at  a 
given  £rsheil  and  energy  and  to  fit  them  with  the  theoretical  model  developed  in 
this  study.  It  also  generates  flux  curves  using  the  fitted  functions  at  predetermined 
times  to  compare  with  the  raw  data  (Ref.  34). 

The  first  two  programs  were  written  with  limited  input  from  the  author.  The  third  was 
written  under  close  and  extensive  coordination  with  this  author,  as  discussed  in  this  chapter. 


The  DIABASE  program,  written  by  Miss  Cheriae  Jarrett  of  Computer  Sciences 
Corporation,  rends  the  raw  data  files,  in  various  formats  created  by  Pfltrer  (Ref.  37),  satel¬ 
lite  by  satellite,  and  stores  them  in  a  three-dimensional  array.  The  points  are  then  ordered 
(within  each  energy  group)  by  increasing  time  since  burst.  To  generate  a  file  of  points  at  a 
specified  /rvalue,  each  data  point  Rvalue  is  compared  to  the  desired  L  and  is  recorded  in 
the  output  file  if  a  match  occurs.  Additional  points  are  ‘created*  by  linear  interpolation  of 
two  successive  data  points  if  they  fall  on  both  sides  of  the  desired  L  and  if  the  data  points 
differ  in  time  by  no  more  than  a  specified  value.  The  interpolation  is  performed  in  L,  B, 
time,  flux,  and  error.  The  output  file  can  be  restricted  to  a  specified  time  ‘window*  and  to 
a  specified  range  of  B  values,  if  desired;  however,  the  data  were  analysed  over  the  full  range 
of  B  and  over  all  times  which  existed  in  the  Data  Base.  The  program  has  tl  e  additional 
capability  to  “0-average"  data  points  at  the  same  L  if  their  individual  B  values  are  within 
a  user  specified  limit  of  each  other  and  they  are  within  the  required  time  difference.  This 
procedure  has  the  effect  of  eliminating  ‘double*  points,  and  was  used  for  the  majority  of 
this  study  with  a  limit  of  one  percent.  Another  capability  of  the  program  is  to  interpolate 
to  a  specified  0-value  in  a  manner  analogous  to  the  //-interpolation  described  above;  this 
capability  was  unused  in  this  study  (Ref.  15:1-2). 


TABLE  D 


Satellites  and  Their  Orbital  Parameters 


Name 

Launch 

Date 

Failure 

Date 

Decay 

Date 

Period 

(min.) 

Perigee 

(mi.) 

Apo-^ee 

(mi.) 

Incl. 

(deg.) 

Akraette  1 

9-28-02 

unk. 

— 

105.4 

620 

634 

80.5 

Alpha 
Upsilon  1 

9-1-82 

— 

10-26-64 

94.4 

189 

4H 

82.8 

Explorer  IS 

10-27-62 

2-9-63 

— 

312.0 

194 

10740 

18.0 

Injun  1 

6-29-61 

3-6-63 

— 

103.8 

534 

631 

67.0 

TVaac 

11-15-01 

7-62 

— 

105.6 

562 

720 

32.4 

Telstar 

7-10-62 

2-21-63 

— 

157.8 

593 

3503 

44.8 

Explorer  4 

7-26-58 

10-6-58 

10-23-59 

110.1 

163 

1372 

50.1 

(Ref.  44:41  45) 
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In  order  to  determine  the  optimnm  time  difference  to  use  Tor  interpolation  acrom  the 
dssirsd  L-vahje,  the  DTABASG  program  wai  modified  under  direction  of  this  antbor.  The 
program  now  presents  a  summary  table  of  the  raw  data  points  for  each  satellite,  showing  the 
frequency  of  oecuranee  of  time  differences  between  successive  data  points  from  Q.l  minute 
to  5.0  minutes,  by  tenths  of  a  minutes.  The  cumulative  percentage  of  points  with  a  specified 
time  difference  or  less  is  also  shown.  These  summaries  are  presented  for  the  bursts  and 
satellites  examined  aa  Tables  III  through  X1D  (Ref.  15.R7-70).  Perusal  of  these  tables  showed 
that  the  optimum  time  difference  for  interpolation  was  2.2  minutes,  primarily  to  include  the 
maximum  number  of  points  in  the  Starfish  Telstar  file.  This  time  difference  was  used  for  ail 
data  examined  in  this  study. 


A  final  capability  of  DTABASG  is  to  add  data  points  to  a  specified  f/- value  file  by 
rounding  points  in  L  to  the  desired  value.  This  procedure  assumes  that  poinu  which  are 
closer  to  the  desired  L  than  some  specified  limit  may  be  considered  to  occur  ai  the  desired 
L.  Within  the  errors  in  satellite  measurements,  this  procedure  is  reasonable,  and  in  fact 
appears  (by  data  comparison)  to  have  been  used  by  Roberts  (Ref.  41:310),  although  he  does 
not  so  explicitly  state.  This  capability  exists  in  the  program,  but  was  not  used  n  this  study 
(Ref.  15). 


The  plotting  program  was  developed  by  Mr.  John  Burgio  of  AFWL  to  generate  all  the 
plots  of  raw  data  and  fitted  curves  shown  in  this  report.  The  program  plots  raw  data  points 
beginning  at  a  specified  start  time  since  burst  (in  days)  and  covering  the  time  period  in  a 
specified  time  “window”  (shown  in  Figure  headings  as  “TW”).  The  L- value,  equatorial  0- 
valne,  and  0-cutoff  values  are  also  shown.  If  the  0-range  of  points  is  restricted,  the  minimum 
and  maximim  values  of  B  are  shown.  The  equatorial  B  value  is  computed  from  L  by  the 

„  0.308 

B°~  i>  ’ 

The  0-cutoff  value  was  initially  computed  by  using  the  dipole  formula  (Ref.  311:55) 


=  J — . . =  VWft-.). 

However,  R  was  later  determined  that  this  approximation  might  be  less  than  optimum;  so 
a  table  of  /rvalues  versus  Bnt  was  used  and  was  linearly  interpolated  across  the  tabulated 
/rvalues  to  find  the  B,nt  value.  This  table  (Table  XIV)  was  prepared  by  AFWL  from  the 
48- term  Jensen-Cain  model  of  the  magnetosphere  (Ref.  IS).  It  was  found  that  use  of  Table 
XIV  produced  Btu ,  (or  yCtt<  or  xeut)  values  more  indicative  of  the  global  cutoff  values.  This 
table  was  used  thereafter  in  both  the  plotting  program  and  the  fitting  program  (discussed 
below).  This  procedure  is  reasonable  in  view  of  the  use  of  the  same  Jensen-Cain  model  to 
produce  the  B~L  coordinates  for  the  Pfitser  Data  Base  (Ref.  37). 


The  plotting  program  plots  flux  as  a  function  of  *  =  y/\  -  (Bo/B).  It  also  has  the 
capability  to  plot  flux  data  points  versus  time  over  all  values  of  L  and  0,  or  flux  versus 
energy  at  specified  0  or  L.  Next  to  the  satellite  name  on  each  plot  which  has  fitted  curves 
is  a  number  in  square  brackets  which  corresponds  to  the  serial  number  of  tbe  run  of  the 
fitting  program  which  prodeuced  the  curves. 


This  author's  input  to  the  plotting  program  was  limited  to  supplying  header  information 
and  specifying  layout  and  format. 

The  third  program  used  in  data  analysis  is  “Program  Rlectrofit,”  written  by  Mr.  Tlarry 
Murphy  of  AFWL  with  the  close  coordination  of  this  author  (Ref.  34).  Electro, It  is  a  least- 
squares  fitting  program  which  uses  the  functions  developed  by  this  author  (Appendix  C) 
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to  (It  data  points  which  are  output  from  DTABASE  in  order  to  compute  the  characteristic 
decay  times  (r„)  and  linear  amplitude  constants  (a*).  These  decay  times  and  constants  are 
determined  at  the  specific  energy  and  Rvalue  which  are  input  with  the  data. 


Table  IV 


$ 


COUNTS  OF  TIME  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHANNELS  OF  SATELL I TE ( S > :  EXPL0RER4 
FOR  BURST:  ARGUS  1.2  .3 


Tl« 

BIN 

NO  DELTA  TIMES 

CUMULATIVE 

(MIN) 

IN  BIN 

PERCENTAGE 

0 

0 

0  1 

442 

63  7 

0 

1 

0  2 

133 

82  9 

0 

2 

0  3 

46 

89  5 

0 

3 

0  4 

10 

90  9 

0 

4 

0  5 

8 

92  1 

0 

5 

0  6 

0 

92  1 

0 

6 

0  7 

1 

92  2 

0 

7 

0  B 

1 

92  4 

0 

8 

0  9 

0 

92  4 

o 

9 

1. 0 

0 

92  4 

1 

0 

1.  1 

0 

92  4 

1 

1 

t  2 

0 

92  4 

1 

2 

1  3 

0 

92  4 

1 

3 

1. 4 

0 

92  4 

1 

4 

1  5 

0 

92  4 

1 

5 

1. 6 

0 

92  4 

1 

6 

1.  7 

0 

92  4 

1 

7 

1.8 

o 

92  4 

1 

8 

1.  9 

0 

92  4 

1 

9 

2.  0 

0 

92  4 

2 

0 

2.  1 

0 

92.  4 

2 

1 

2  2 

0 

92.  4 

2 

2 

2.  3 

0 

92  4 

2 

3 

2.  4 

1 

92  5 

2 

4 

2  5 

0 

92  5 

2 

5 

2  6 

0 

92  5 

2 

6 

2  7 

0 

92.  5 

2 

7 

2  8 

1 

92  7 

2 

a 

2.  9 

•t 

92  8 

2 

9 

3  0 

0 

92  8 

3 

0 

3.  1 

1 

92  9 

3 

1 

3  2 

1 

93  1 

3 

2 

3  3 

0 

93  | 

3 

3 

3.  4 

0 

93  1 

3 

4 

3  5 

0 

93  1 

3 

5 

3.  6 

0 

93  1 

3 

6 

3  7 

0 

93  1 

3 

7 

3  8 

0 

93.  1 

3 

8 

3  9 

0 

93  1 

3 

9 

4.  0 

0 

93  1 

0 

4  1 

0 

93  1 

1 

4.  2 

0 

93  1 

2 

4.  3 

0 

93  1 

3 

4.  4 

0 

93  1 

4 

4  S 

0 

93  1 

5 

4  & 

0 

93.  1 

.  6 

4  7 

0 

93  1 

7 

4  8 

0 

93  1 

8 

4.  9 

0 

93  1 

.  9 

5  0 

0 

93  1 

S. 

o 

A 

R 

48 

too  0 

Table  V 


3 

jj  nr 


COUNTS  OF  TIME  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHAhWELS  OF  SATELL I TE < S ) :  ALLOUETTE1 
FOR  BURST  RUSSIAN  1 


TI* 

BIN 

NO.  DELTA  TIMES 

CUMULATIVE 

(MIN> 

IN  BIN 

PERCENTAGE 

0  0 

0.  1 

0 

0  0 

0.  1 

0.  2 

3413 

72  7 

0  2 

0.  3 

0 

72  7 

0.  3 

0  4 

972 

84  8 

0.  4 

0.  9 

180 

88  7 

0  S 

0.  6 

15 

89  0 

0  6 

0.  7 

105 

91  2 

0.  7 

0.  8 

0 

91  2 

0  8 

0  9 

46 

92.  2 

0  9 

1 . 0 

46 

93  2 

1  0 

1  1 

1 

93  2 

1.  1 

1.  2 

5 

93  3 

1  2 

1.  3 

0 

93  3 

1. 3 

1.  4 

15 

93  6 

1. 4 

1. 9 

4 

93  7 

1.  5 

1. 6 

3 

93  8 

1.6 

1.  7 

5 

93  9 

1.7 

1.8 

0 

93.  9 

1  8 

1.  9 

12 

94  1 

1  9 

2.  0 

1 

94.  2 

2.  0 

2.  I 

1 

94  2 

2  1 

2.  2 

2 

94  2 

2.  2 

2.  3 

0 

94  2 

2  3 

2.  4 

2 

94.  3 

2.  4 

2.  9 

1 

94.  3 

2  9 

2.  6 

0 

94  3 

2.  6 

2.  7 

3 

94  4 

2  7 

2  8 

0 

94  4 

2  8 

2  9 

0 

94  4 

2  9 

3.  0 

0 

94  4 

3  0 

3.  1 

1 

94.  4 

3  1 

3  2 

0 

94  4 

3  2 

3.  3 

0 

94.  4 

3  3 

3.  4 

3 

94  4 

3  4 

3  9 

0 

94  4 

3  9 

3.  6 

0 

94  4 

3  6 

3  7 

0 

94  4 

3.  7 

3  8 

0 

94  4 

3.  8 

3  9 

1 

94  5 

3  9 

4  0 

1 

94  5 

4  0 

4.  1 

0 

94  5 

4  1 

4.  2 

3 

94  5 

4  2 

4.  3 

0 

94  5 

4.  3 

4.  4 

1 

94  6 

4.  4 

4.  9 

0 

94  6 

4.  9 

4  6 

0 

94  6 

4.  6 

4  7 

0 

94  6 

4  7 

4.  8 

0 

94  6 

4.  8 

4  9 

0 

94  6 

4.  9 

9.  0 

0 

94  6 

9. 

O 

A 

■ 

259 

100  0 
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Table  VI 


M  re 


COUNTS  OP  TIME  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHANNELS  OF  SATELLITES):  EXPLORER  15 

FOR  BURST  RUSSIAN  1 

T IPE  BIN  NO.  DELTA  TIMES  CUMULATIVE 

(MIN)  IN  BIN  PERCENTAGE 


0  0  0  1 

0.  1  0  2 

0.2  0  3 

0.  3  0  4 

0  4  0  5 

0.5  0  6 

0  6  0  7 

0.7  OB 
OB  0  9 
O.  9  10 

10  11 
1.1  12 

12  13 

13  1.4 

1.4  1.9 

15  16 

1.6  1.7 

1.7  18 

1.  B  1.9 

1.9  2  0 

2  0  2  1 

2.1  2.2 
2  2  2  3 

2  3  2  4 

2  4  2  5 

2  5  2  6 

2.  6  2  7 

2  7  2  8 

2  8  2  9 

2.9  3.0 

3.  0  3  1 

3  1  3  2 

3  2  3  3 

3  3  3  4 

3  4  3  5 

3.  5  3  6 

3  6  3  7 

3.7  3  8 

3  8  3  9 

3  9  4  0 

4.0  4  1 

4  1  4  2 

4  2  4  3 

4.3  4  4 

4.4  4  5 

4.5  4  6 

4  6  4.7 

4  7  4  8 

4  8  4  9 

4  9  5  0 

5  0<« 


0 

0 

O 

2 

12 

0 

5 

6 
1 
O 

13 

8 

0 

1 

O 

0 

11 

3 

O 

o 

0 

3 

4 
O 
0 
o 
0 
0 
1 
0 
2 
0 
0 
3 
7 
0 
3 
0 
0 
1 
0 
0 
0 
o 
0 
1 
o 
0 
0 
0 
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0  0 
0  0 
0  0 
1  7 
12  0 
12  0 
16  2 

21  4 

22  2 
22  2 
33.  3 
40  2 

40  2 

41  O 
41  O 
41  0 
50  4 
53  0 
53  0 
53  0 
53  0 
55  6 
59  O 
59  0 
59  0 
59.  0 
59.  0 
59  0 
59  B 
59  8 
61  5 
61  5 
61.  5 
64.  1 
70  1 
70  1 
72  6 
72  6 

72  6 

73  5 
73  5 
73  5 
73  5 
73  5 

73  5 

74  4 
74  4 
74  4 
74  4 
74  4 

100  0 
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Table  VTI 


COUNTS  or  TIME  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHANNELS  OF  SATELL I TE < S > :  TELSTAR 
FOR  BURST  RUSSIAN  1 


TIME  BIN 
(MIN) 


0  0 

0  1 

0  1 

0  2 

0  2 

0  3 

0  3 

0  4 

0  4 

0  5 

0  5 

0  6 

0  6 

0  7 

0  7 

0  8 

0  8 

0.  9 

0  9 

1  O 

1  0 

1.  1 

1  1 

1  2 

1  2 

1  3 

1  3 

1  4 

1. 4 

1.  5 

1  5 

1. 6 

1. 6 

1  7 

1  7 

1. 8 

1  8 

1  9 

1  9 

2  0 

2  O 

2.  1 

2  1 

2  2 

2.  2 

2.  3 

2.  3 

2  4 

2.  4 

2  5 

2  9 

2  6 

2.  6 

2  7 

2.  7 

2  8 

2.  8 

2.  9 

2  9 

3  0 

3.  0 

3  1 

3.  1 

3.  2 

3  2 

3  3 

3  3 

3  4 

3  4 

3  5 

3  9 

3.  6 

3.  6 

3  7 

3  7 

3  B 

3  8 

3  9 

3  9 

4  0 

4.  0 

4  1 

4.  1 

4  2 

4.  2 

4  3 

4  3 

4  4 

4  4 

4  9 

4  9 

4  6 

4  6 

4  7 

4  7 

4.  8 

4.  8 

4  9 

4.  9 

9  0 

9. 

o 

A 

1 

DELTA  TIMES 
IN  BIN 

0 

0 

0 

0 

1 

O 

0 

0 

0 

7 

1 

0 

0 

0 

o 

0 

o 

1 

2 

3637 

3 

3 

0 

0 

0 

0 

0 

0 

0 

6 

0 

0 

0 

0 

0 

0 

0 

0 

1 

134 

0 

0 

O 

2 
0 
0 
0 
0 
0 
9 

413 


CUMULATIVE 

PERCENTAGE 

O  0 
0  0 
O  0 
O.  0 
O  0 
0  0 
0  O 
0  0 
0  O 
0  2 
O.  2 
0  2 
0  2 
0  2 
0  2 
0  2 
0  2 
0  2 
O  3 
86.  6 
86  6 
86  7 
86  7 
86  7 
86  7 
86  7 
86  7 
86  7 
86  7 
86.  8 
86  8 
86  8 
86  8 
86  8 
86  8 
86  8 
86  8 
86  8 
86  9 
90  0 
90  0 
90  0 
90  0 
90  1 
90  1 
90  1 
90  1 
90  t 
90  1 
90  2 
100  0 


Table  VIII 


.1 


9  r* 


■ 


COUNTS  OF  TIME  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHANNELS  OF  SATELL I TE < S ) :  ALLOUETTE1 
FOR  BURST  RUSSIAN  2 


TIME  BIN 
(MINI 


NO.  DELTA  TirCS 
IN  BIN 


CUMULATIVE 

PERCENTAGE 


l 


i 


3 


1 


■ 


« 


0.  0 

0 

i 

45 

1 

1 

0.  1 

0 

2 

3107 

78 

9 

0  2 

0 

3 

0 

78 

o 

0.  3 

0 

4 

409 

89 

1 

0  4 

0 

5 

91 

91 

4 

0  3 

0 

A 

7 

91 

5 

0  6 

0 

7 

32 

92 

8 

0  7 

0 

8 

0 

92 

0 

0  8 

o 

9 

29 

93 

6 

0  9 

I 

0 

25 

94 

2 

1  0 

1 

1 

3 

94 

3 

1.  1 

1 

2 

13 

94 

A 

1  2 

1 

3 

0 

94 

A 

1.  3 

1 

4 

7 

94 

8 

1. 4 

I 

5 

1 

94 

8 

1.  3 

1 

A 

0 

94. 

8 

1. 6 

1 

7 

A 

94 

9 

1.  7 

I 

8 

0 

94 

9 

1  8 

1 

9 

4 

95 

0 

I.  9 

2 

0 

1 

95 

1 

2  0 

2 

1 

0 

95 

1 

2  1 

2 

2 

4 

95 

2 

2  2 

2 

3 

0 

95 

2 

2.  3 

2 

4 

1 

95 

2 

2  4 

2 

5 

3 

95 

3 

2  3 

2 

A 

0 

95 

3 

2  & 

2 

7 

0 

95 

3 

2  7 

2 

8 

0 

95 

3 

2  8 

2 

9 

0 

95. 

3 

2  9 

3 

0 

0 

95 

3 

3.  0 

3 

1 

0 

95 

3 

3  1 

3 

2 

0 

95 

3 

3.  2 

3 

3 

0 

95 

3 

3.  3 

3 

4 

0 

95 

3 

3.  4 

3 

5 

0 

95 

3 

3  3 

3 

A 

0 

95 

3 

3  6 

3 

7 

0 

95 

3 

3  7 

3 

8 

0 

95 

3 

3.  8 

3 

9 

0 

95 

3 

3  9 

4 

0 

0 

95 

3 

4  0 

4 

I 

0 

95 

3 

4  1 

4 

2 

0 

95 

3 

4.  2 

4 

3 

0 

93 

3 

4  3 

4 

.  4 

0 

95 

3 

4  4 

4 

3 

0 

95 

3 

4  3 

4 

A 

0 

93 

3 

4.  A 

4 

7 

0 

93 

3 

4.  7 

4 

8 

0 

93 

3 

4  8 

4 

9 

0 

95 

3 

4  9 

3 

0 

0 

95 

3 

3 

o 

A 

a 

189 

100 

0 

Table  DC 


COUNTS  OF  Tire  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHANNELS  OF  SATELL I TE ( S ) :  EXPLORER  15 

FOR  BURST:  RUSSIAN  2 


Tire 

BIN 

NO  DELTA  TIMES 

CUMULATIVE 

(MIN) 

IN  BIN 

PERCENTAGE 

0  0 

0.  1 

0 

0  0 

0.  1 

0.  2 

0 

0  0 

0  2 

0.  3 

0 

O  0 

0.  3 

0.  4 

71 

2  9 

0.  4 

0  5 

264 

13  8 

0  5 

0.  6 

344 

28  0 

0  6 

0  7 

142 

33  9 

O  7 

0.  8 

77 

37  0 

0  B 

0  9 

31 

38  3 

0  9 

1.0 

83 

41  7 

1  0 

1.  1 

65 

44  4 

i.  1 

1.2 

68 

47  2 

12 

1. 3 

64 

49  9 

1  3 

1.  4 

63 

52  5 

1.  4 

1.  9 

115 

57  2 

1  5 

1. 6 

75 

60  3 

1.6 

1.  7 

129 

65  4 

1  7 

1.  8 

36 

66  9 

1  e 

1  9 

11 

67  4 

19 

2.  0 

33 

68  7 

2  0 

2  1 

12 

69  2 

2  1 

2  2 

29 

70  4 

2  2 

2  3 

32 

71  8 

2  3 

2  4 

17 

72  5 

2  4 

2.  9 

42 

74  2 

2  9 

2  6 

14 

74  8 

2  6 

2  7 

26 

75  8 

2  7 

2  8 

11 

76  3 

2  8 

2  9 

16 

76  9 

2  9 

3.  0 

10 

77  4 

3  0 

3  1 

10 

77  8 

3.  1 

3  2 

10 

78  2 

3  2 

3  3 

7 

78  5 

3  3 

3  4 

8 

78  8 

3  4 

3  9 

22 

79  7 

3  9 

3  6 

23 

80  7 

3  6 

3  7 

45 

82  5 

3  7 

3  8 

14 

83  1 

3.  8 

3.  9 

18 

83  8 

3  9 

4.  0 

27 

84  9 

4  0 

4  1 

15 

85  6 

4  1 

4  2 

10 

86.  0 

4.  2 

4.  3 

6 

86  2 

4  3 

4.  4 

10 

86  6 

4  4 

4  9 

4 

86  8 

4.  9 

4.  6 

19 

87  4 

4  6 

4  7 

10 

87  8 

4.  7 

4  8 

4 

88  0 

4  8 

4  9 

14 

88  6 

4  9 

9  0 

8 

88  9 

9 

o 

A 

t 

269 

100  0 

.V, 


Table  X 


COUNTS  OF  TINE  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHANNELS  OF  SATELL ITE < S ) :  TEL STAR 
FOR  BURST:  RUSSIAN  2 

TirC  BIN  NO.  DELTA  TlfCS  CUMULATIVE 

(NIN)  IN  BIN  PERCENTAGE 


0  0 

0  1 

0 

0  0 

0.  1 

0  2 

0 

0  0 

0.  2 

0.  3 

0 

0  0 

0.  3 

0.  4 

0 

0  0 

0.  4 

0  5 

0 

0  0 

0  5 

0  6 

0 

0  0 

0  6 

0  7 

0 

0  0 

0  7 

0  8 

0 

0  0 

0  8 

0  9 

0 

0.  o 

0  9 

1. 0 

2 

0  1 

1.0 

1.  1 

0 

0  1 

1.  1 

1. 2 

0 

0  1 

12 

1. 3 

0 

0  1 

1. 3 

1. 4 

0 

0  1 

1.  4 

1.  5 

0 

0  1 

1.  9 

16 

0 

0  1 

1.  A 

1. 7 

0 

0  1 

1  7 

1.8 

1 

0  1 

1  B 

1.  9 

0 

0.  1 

1.9 

2.  0 

3244 

83.  6 

2  0 

2.  1 

49 

87.  3 

2  1 

2.  2 

1 

87  3 

2  2 

2  3 

0 

87  3 

2.  3 

2  4 

0 

87  3 

2.  4 

2.  5 

0 

87  5 

2.  3 

2.  6 

0 

87  3 

2.  6 

2.  7 

0 

87.  3 

2.  7 

2.  8 

0 

87  3 

2  8 

2  9 

0 

87.  3 

2  9 

3  0 

3 

87  6 

3  0 

3  1 

0 

87  6 

3.  1 

3  2 

0 

87.  6 

3  2 

3  3 

0 

87  6 

3  3 

3  4 

0 

87  6 

3.  4 

3  3 

0 

87  6 

3.  3 

3  6 

0 

87.  6 

3  6 

3  7 

0 

87  6 

3  7 

3  8 

0 

87  6 

3.  8 

3.  9 

2 

87  7 

3.  9 

4  0 

100 

91  5 

4.  0 

4.  1 

2 

91  6 

4  1 

4  2 

0 

91  6 

4  2 

4  3 

0 

91  6 

4  3 

4.  4 

0 

91  6 

4.  4 

4  3 

0 

91  6 

4  3 

4  6 

0 

91.  6 

4.  6 

4.  7 

0 

91  6 

4  7 

4  8 

0 

91  6 

4  0 

4  9 

0 

91  6 

4  9 

9  0 

0 

91  6 

3. 

■ 

V 

o 

221 

100  0 

no 


Tabic  XI 


nr 


COUNTS  OF  TINE  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHANNELS  OF  SATELL I TE ( S ) :  ALLOUETTE1 
FOR  BURST:  RUSSIAN  3 


TI*C 

BIN 

NO  DELTA  TIMES 

CUMULATIVE 

(MIN) 

IN  BIN 

PERCENTAGE 

0  0 

0  1 

123 

0  8 

0  1 

0  2 

10901 

74  0 

0  2 

0.  3 

49 

74  3 

0  3 

0.  4 

1995 

85  0 

0  4 

0.  9 

491 

88  3 

0  9 

0.  6 

95 

88  7 

0  6 

0  7 

289 

90  6 

0  7 

0.  8 

11 

90  7 

0  B 

0  9 

128 

91  5 

0  9 

1. 0 

74 

92  0 

10 

1  1 

8 

92  1 

1  1 

1  2 

95 

92  5 

1.2 

1.  3 

3 

92  5 

1. 3 

1.  4 

62 

92  9 

1.4 

1.  9 

28 

93  1 

1  5 

1.  & 

9 

93  1 

16 

1.  7 

18 

93  3 

1.7 

1.  8 

0 

93  3 

18 

1.9 

18 

93  4 

1.9 

2.  0 

15 

93  5 

2.  0 

2.  1 

2 

93  5 

2.  1 

2.  2 

10 

93  6 

2  2 

2.  3 

0 

93  6 

2  3 

2.  4 

13 

93  7 

2.  4 

2.  9 

7 

93  7 

2.  9 

2.  6 

1 

93  7 

2.  6 

2.  7 

7 

93  8 

2  7 

2  8 

0 

93  8 

2  B 

2.  9 

'4 

93  8 

2.  9 

3.  0 

2 

93  8 

3  0 

3.  1 

2 

93  8 

3.  1 

3  2 

9 

93  8 

3  2 

3  3 

0 

93  8 

3.  3 

3.  4 

2 

93  9 

3  4 

3.  9 

0 

93  9 

3  9 

3.  6 

1 

93  9 

3.  6 

3.  7 

1 

93  9 

3  7 

3  8 

0 

93  9 

3.  8 

3  9 

1 

93.  9 

3.  9 

4.  0 

0 

93  9 

4.  0 

4.  1 

0 

93  9 

4  1 

4  2 

0 

93.  9 

4.2 

4  3 

0 

93  9 

4.  3 

4.  4 

0 

93  9 

4  4 

4.  9 

0 

93  9 

4  9 

4.  6 

1 

93  9 

4.  6 

4.  7 

0 

93  9 

4  7 

4.  8 

0 

93  9 

4  8 

4.  9 

'  0 

93  9 

4.  9 

9.  0 

0 

93  9 

9 

O 

A 

■ 

911 

100  0 

Table  XII 


COUNTS  OF  TINE  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHANNELS  OF  SATELL I TE  ( S > :  EXPLORER  15 
FOR  BURST:  RUSSIAN  3 

TINE  BIN  NO.  DELTA  TICES  CUMULATIVE 

<NIN)  IN  BIN  PERCENTAGE 


0  0 

0  1 

46 

0.  3 

0.  1 

0  2 

0 

0  3 

0.  2 

0  3 

77 

0  8 

0.  3 

0.  4 

988 

4  8 

0.  4 

0  5 

1654 

16.  0 

0.  5 

0.  6 

2144 

30.  5 

0.  6 

0.  7 

1629 

41.  5 

0.  7 

0.  8 

970 

48  0 

o  e 

0.  9 

274 

49.  9 

0  9 

1.  0. 

194 

51.  2 

10 

1.  1 

226 

52  7 

1.  1 

1.  2 

210 

54  1 

1.2 

1.  3 

293 

56.  1 

1.3 

1.4 

273 

58  0 

1.4 

1.  9 

474 

61  2 

1  5 

1.  6 

442 

64  1 

1.6 

1.  7 

618 

68  3 

17 

1.  8 

313 

70  4 

1.8 

1.  9 

194 

71.  7 

1.9 

2  O 

131 

72.  6 

2.  0 

2.  1 

B2 

73  2 

2.  1 

2.  2 

92 

73  8 

2.  2 

2.  3 

74 

74  3 

2.  3 

2.  4 

94 

74  7 

2.  4 

2.  9 

45 

75  0 

2.  9 

2  6 

35 

75  2 

2.  6 

2.  7 

93 

75  6 

2.  7 

2.  8 

34 

75  8 

2.  8 

2.  9 

37 

76  0 

2  9 

3.  0 

40 

76  3 

3.  0 

3.  1 

29 

76  5 

3  1 

3.  2 

47 

76.  8 

3.  2 

3  3 

52 

77  2 

3.  3 

3.  4 

66 

77  6 

3.  4 

3.  5 

93 

78  0 

3  9 

3.  6 

79 

78  5 

3.  6 

3  7 

195 

79  6 

3  7 

3.  8 

91 

80.  2 

3.  S 

3.  9 

127 

81.  0 

3.  9 

4.  0 

111 

81  8 

4.  0 

4.  1 

83 

82  3 

4.  1 

4  2 

117 

83  1 

4  2 

4  3 

87 

83.  7 

4.  3 

4.  4 

110 

84.  5 

4.  4 

4  9 

73 

89  0 

4  9 

4.  6 

100 

85  6 

4.  6 

4.  7 

76 

86  2 

4.  7 

4  8 

63 

86  6 

4  8 

4.  9 

80 

87  1 

4  9 

9.  0 

59 

87  5 

9. 

0<- 

1848 

100  0 

.*** 


R  « 


Table  XID 


COUNTS  OF  TIME  DIFFERENCES  BETWEEN  NONZERO  RAW  DATA  POINTS 
FOR  ALL  CHAMMELS  OF  SATELLITES):  TELSTAR 
FOR  BURST:  RUSSIAN  3 


TlfC  BIN 
(HIN) 


NO. 


DELTA  TDCS 
IN  BIN 


CUMULATIVE 
PER  CENT  ACE 


0.  0 

0.  1 

3 

0  0 

0  1 

0  2 

0 

0  0 

0  2 

0  3 

0 

0  0 

0.  3 

0.  4 

0 

0  0 

0  4 

0  5 

0 

0  0 

0  5 

0  A 

0 

0  0 

0.  6 

0.  7 

0 

0  0 

0.  7 

0  8 

0 

0  0 

0  B 

0  9 

1 

0  0 

0:  9 

1.  0 

64 

0  3 

10 

1.  1 

4 

0  3 

1. 1 

1.2 

2 

0  3 

1.2 

1.  3 

0 

0  3 

1  3 

1  4 

0 

0  3 

1. 4 

1. 3 

1 

0  3 

1.  S 

16 

3 

0  3 

1.  6 

1.  7 

3 

0  3 

1.  7 

1.8 

6 

O  4 

1.  B 

1.9 

15 

0  4 

1  9 

2.  0 

20107 

82  9 

2.  0 

2.  1 

121 

83  4 

2.  1 

2.  2 

7 

83  3 

2  2 

2  3 

6 

83  5 

2  3 

2  4 

1 

83  5 

2  4 

2.  3 

0 

83  3 

2  3 

2  6 

0 

83  3 

2  & 

2  7 

0 

83  5 

2  7 

2.  8 

0 

83  5 

2.  B 

2  9 

0 

83  5 

2.  9 

3  0 

25 

83  6 

3  0 

3  1 

1 

83  6 

3.  1 

3.  2 

3 

83  6 

3  2 

3  3 

1 

83  6 

3.  3 

3.  4 

0 

83  6 

3.  4 

3.  3 

0 

83  6 

3  3 

3.  & 

0 

83  6 

3.  6 

3.  7 

3 

83  6 

3.  7 

3.  8 

2 

83  6 

3.  8 

3  9 

9 

83  7 

3.  9 

4  0 

1027 

87  9 

4.  0 

4  1 

5 

87  9 

4.  1 

4.  2 

1 

87  9 

4.  2 

4.  3 

0 

87  9 

4  3 

4  4 

3 

87  9 

4.  4 

4  3 

1 

87  9 

4  3 

4  6 

0 

87  9 

4.  6 

4.  7 

0 

87  9 

4.  7 

4  8 

0 

87  9 

4.  8 

4.  9 

1 

87  9 

4.  9 

3.  0 

12 

88  0 

3. 

OO 

2929 

100  0 

ns 


Table  XTV 

Global  Magnetic  Cutoff  Values 


Value 

1.0949 

0.23901 

1.1152 

.22998 

1.1356 

.22482 

1.1550 

.22288 

1.1782 

.22353 

1.2628 

,240t7 

1.4817 

.26770 

1.8331 

.27573 

2.3168 

.29762 

2.9330 

.33311 

3.6816 

.37001 

4.5626 

.40274 

5.5760 

.43178 

6.7218 

.45502 

8.0000 

.46951 

Note:  Tbe  Table  may  be  interpolated  for  other  L-rahns  (Ref.  19). 


r? 


The  overall  flow  for  the  program  is  as  follows  (see  Appendix  E): 

(1)  Use  an  initial  guess  for  the  fundamental  decay  time  and  the  theoretical  ratios 
of  the  higher-mode  de--y  times  to  compute  the  starting  decay  tin  es.  The  first 
guess  at  fundamental  decay  time  is  shown  in  Table  XV  as  a  function  of  L  and 
energy.  This  table  is  interpolated  in  L  and  R  for  starling  value. 

(2)  Numerically  compute  the  integrals  of  the  functions  shown  in  Appendix  C  at  each 

data  point  using  the  initial  guessed  decay  times.  The  integration 's  performed 
by  a  12-point  Oauss/Legcndre  Quadrature  at  each  data  point  for  each  function. 
Only  the  fundamental  decay  mode  and  two  higher  modes  were  cc  mputed. 

(3)  Compute  the  linear  amplitude  constants  using  a  least-squares  matrix-solving 
routine.  Compute  the  standard  error  of  the  fit.  The  individiml  points  are 
weighted  by  the  reciprocal  of  their  fractional  errors,  which  have  been  normal¬ 
ised  to  an  average  value  of  one. 

(4)  Optimise  the  initially  computed  j/eut  value  by  using  a  “Golden  Section  Min- 
imim”  function.  This  section  attempts  to  let  the  data  determine  the  “best" 
cutoff  value.  If  the  data  are  not  distributed  down  toward  the  loss  cone,  tbe 
Colden  Minimum  function  will  provide  an  erroneous  ye„t  value.  Hence,  the 
function  specifies  that  the  computed  yeut  must  be  within  ±10  percent  of  the 
initial  input  value.  If,  in  the  iteration  process,  tbe  function  tries  to  exrccd  these 
limits,  the  initial  guess  of  yeut  is  specified  as  the  correct  value. 

(5)  Repeat  the  steps  (2)  and  (3)  with  the  optimum  y,mt  value 

(8)  To  determine  the  best  combination  of  decay  times  and  amplitude  constants, 
the  decay  times  arc  changed  by  a  “random-walk"  method,  either  increased 
or  decreased.  The  amount  by  which  the  decay  times  may  randomly  vary  is 
determined  by  using  a  cumulative  normal  probability  distribution  about  the 
initial  value. 

The  standard  deviation  used  was  estimated  according  to  this  author’s  best  guess  of 
the  accuracy  of  the  tables  and  curves  shown  by  Stassinopoulos  (Ref.  53:31-32,  40-44)  and 
West  (Ref.  83:50-54).  Tbe  Stassinopoulos  tables  do  not  always  correlate  with  each  other 
or  with  the  Stassinopoulos  curves,  and  the  West  curves  do  not  always  correlate  with  the 
Stassinopoulos  in  the  ^-regions  of  overlap.  The  initial  fundamental  decay  times  in  Table  XV 
were  computed  by  this  author  by  interpolation  of  the  tables  and  curves  in  Reference  53  and 
by  reading  and  extrapolation  of  the  curves  in  Reference  83.  The  author’s  estimates  of  their 
standard  deviations  are  shown  in  Table  XVI.  Of  particular  note  is  the  very  low  confidence 
in  decay  limes  at  high  L  values  (above  2.2)  and  at  high  energies  (above  2.0  MeV).  There 
have  been  few  computations  of  decay  times  in  these  regions,  and  the  results  vary  widely 
(Refs.  53,  53). 

The  effect  of  using  the  cumulative  normal  probability  distribution  and  the  above  stan¬ 
dard  deviations  is  to  restrict  the  “random-walk"  of  the  decay  times  away  from  their  initial 
values.  The  farther  away  from  the  initial  “mean  value"  that  the  random  val  «e  generator 
places  the  new  decay  time,  the  greater  the  "push"  that  is  created  back  toward  he  mean  for 
the  next  walk.  The  distance  away,  however,  is  defined  in  terras  of  the  standaid  deviation, 
so  that  r  values  in  which  there  is  low  confidence  may  vary  by  larger  amounts  than  may  r 
values  in  which  them  is  greater  confidence.  The  r’s  are  also  restricted  by  the  ro  islraint  that 
ro  >  r,  >  r2. 


•U 


rs 

With  each  change  of  decay  times,  steps  (2)  and  (3)  are  repeated  np  to  96  times,  in  an 
attempt  to  get  a  “better*  fit  by  reducing  the  standard  error  of  the  fit. 
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((Tom (lull'd  from  IM  f»3  31  32,  10  It  ;u.<l  K.f.  I»3.  50  51) 


Table  XVT 

Percent  Standard  Deviation  of  Initial 
Fundamental  Decay  Times  (in  Table  XV) 


L-Rango 

Energy  Range 
MeV 

1.02- 

1.20 

1.21- 

1.40 

1.41- 

1.80 

1.81- 

2.20 

> 

2.20 

E  <  0.8 

50 

30 

20 

35 

200 

0.8  <  E  <  2.0 

30 

20 

20 

75 

100 

E  >  2.0 

50 

20 

50 

100 

100 

(estimated  from  Ref.  53:31-32,  40-44  and  Ref.  83:50-54) 


(7)  If  a  better  fit  is  found  in  less  than  98  random  walks,  and  if  less  than  .12  ‘better’ 
▼aloes  hare  been  found,  steps  (2),  (3)  and  (0)  are  repeated  up  to  32  times,  using 
the  latest  “best’  decay  times. 

(8)  Up  to  fire  repeats  of  steps  (4)  and  (5)  may  be  carried  out,  but  only  if  32  “better" 
fits  are  found.  Otherwise,  only  the  first  yeut  optimisation  is  performed. 

(9)  The  program  will  always  end  at  the  point  where  96  random- walk  *  have  failed 
to  reduce  the  error  of  the  fit  from  the  previous  “best-fit’.  The  lecay  times, 
yeut,  and  amplitude  constants  are  printed  for  each  “better"  fit,  v  it h  the  final 
values  being  those  of  the  “best"  fit. 

(10)  After  the  best  fit  is  found,  the  program  uses  the  fitted  constants  and  decay 
times  to  compute  fluxes,  at  specified  times,  from  the  equator  to  the  best  ruled 
value.  These  specified  times  correspond  to  default  times  plus  1/2  ime  window 
in  the  program  which  plots  raw  data.  Thus,  the  computed  fluxe:  are  used  to 
generate  curves  which  are  plotted  over  raw  data  points  at  approximately  the 
same  times. 

Variations  in  the  program  were  tried  in  order  to  study  the  effects  of  the  filling  process. 
The  principal  modifications  were  in  the  random  decay-time  variations. 

The  first  modification  was  to  fix  the  ratios  of  the  r’s  to  their  theoretical  values  of 
1:1/9:1/25  (Ref.  48:16-19).  The  effect  of  this  modification  is  to  significantly  reduce  the 
degrees  of  freedom  available  to  fit  the  data. 

The  second  modification  was  a  compromise  between  the  above  and  complete  freedom  to 
random  walk.  This  variation  allowed  r0  and  r,  to  vary  freely  (subject  to  the  usual  contraint 
that  r0  >  n),  but  rt  and  r2  were  fixed  at  their  theoretical  ratio  of  25:9.  The  principal  reason 
for  this  attempt  was  the  very  short  time-span  covered  by  some  satellite/burst  combinations. 
For  example,  the  Russian  2  Burst  covered  only  4  1/2  days  prior  to  the  Russian  3  injection. 
If  the  data  cover  a  time-span  which  is  small  compared  to  the  fundamental  decay  constant, 
then  the  fitted  value  of  that  decay  time  may  be  suspect.  However,  if  the  time-s  ran  is  not  so 
small  compared  to  the  decay  times  of  the  higher  modes,  one  should  have  mote  confidence 
in  the  fitted  values  of  those  decay  time.  Stated  another  way,  the  fit  may  not  b<-  sensitive  to 
the  fundamental,  but  may  still  be  sensitive  to  the  higher  modes. 

In  this  analysis,  not  all  /rvalues  covered  by  the  satellite  in  Table  II  were  examined. 
Over  1000  plots  of  Dux  versus  z  were  examined  at  various  //-values  for  the  ss.tellite/bursl 
combinations  of  Tables  I  and  II.  In  many  cases,  the  satellite  orbits  were  <>f  such  high 
inclination  (toward  polar)  that  data  were  not  available  close  to  equatorial  ^-values.  Because 
of  this,  and  because  of  the  very  large  volume  of  data,  those  satellite/burst///  combinations 
which  provided  the  most  consistent  coverage  from  the  equator  to  the  loss  cone  were  chosen 
for  study. 

(Generally,  the  data  fell  into  three  regions  of  x-space  (with  some  overlap),  legion  I  was 
from  the  equator  to  the  point  where  llux  dropped  off  sharply  into  the  loss  cone  (i  =  0  to 
z  0.8  or  0.8).  Region  II  was  the  area  where  the  Qux  “turned  the  corner"  from  a  rather 
“horiiontal"  curve  to  a  rather  “vertical"  curve  into  the  loss  cone  (z  ss  0.7  to  0.’.>).  Region  ;d 
was  the  loss  cone  region,  where  the  flux  levels  dropped  sharply  to  the  cutoff  v  lue  ( z  .9 
for  a  typical  t-value). 


Representative  plots  of  raw  data  in  the  three  regions  are  shown  in  Figures  4  through  ft. 
These  plots  are  typical  in  that  data  seldom  cover  all  three  regions  in  a  short  time-span.  Gaps 
in  the  data  appear  because  of  satellite  orbital  coverage  and  because  of  lack  of  on-satellite 
recorders.  When  a  satellite  lacked  a  recorder,  data  were  collected  on '.7  as  the  satellite  passed 
above  the  radio  horizon  of  a  ground  tracking  station.  Figure  'J  is  an  eight-hour  plot  of  typical 
satellite  coverage  over  all  B-L  space. 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM2)/SEC 


BURST  :  RUSSIAN  2 
TIME  :  0  (dy)  0  (hr)  ( 
L  :  2.40  Beq  :  .022  BC 
Bmin  (Gauss)  :  0.001 

- , - r~-  ■  i - 


io5 


io4 


10s  I - * - * - * - ‘ - 

0 

X  =  (1- 

Fifture  4(a).  Raw  Allouette  Data  Collected  Oi 


io9 

10" 

io7 

io6 
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INTEGRAL  ELECTRON  FLUX  (ELECTR0NS/CM2)/SEC 


BURST  :  RUSSIAN  2  SAT  :  ALLOUETTE1 
TIME  :  3  (dy)  0  (hr)  0.00  (min)  3.000(days) 
L  :  2.50  Beq  :  .020  BCUT  :  0.515  TW  (days)  :  1.0 
Bmin  (Gauss)  :  0.001  B max  (Gauss)  :  1.000 


LEGEND 

a  =  3.90  MEV 
o  =  3.90  MEV 
a  =  3.90  MEV 


X  =  (1-Beq/B) ' 

Figure  4(b).  Raw  Allonette  Data  Collected  Only  in  Region  ID  (3  days  posi -burst) 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM2)/SEC 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM^/SEC 


BURST  :  STARFISH  SAT  :  TELSTAR 
TIME  :  15  (dy)  0  (hr)  0.00  (min)  I5.000(days) 
L  :  1.70  Beq  :  .063  BCUT  :  0.460  TW  (days)  :  5.0 
Bmin  (Gauss)  :  0.001  Bmax  (Gauss)  :  1.000 


LEGEND 

=  .220  MEV 
=  .320  MEV 
=  .420  MEV 
=  .660  MEV 


X  =  (1-Beq/B) ' 

Figure  5.  Raw  Telstar  Data  Collected  in  Region  III  (mid-range  time) 


INTECRAL  ELECTRON  FLUX  (ELECTRONS/CVT)/SEC 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM2)/SEC 


BURST  :  STARFISH  SAT  :  TELSTAR 
TIME  :  100  (dy)  0  (hr)  0.00  (min)  lOO.OCO(days) 
L  :  2.00  Beq  :  .038  BCUT  :  0.487  TW  (days)  :  5.0 


0  0.5  1 

X  =  (1-Beq/B)1/2 


Figure  8(b).  Raw  Telstar  Data  in  Region*  I,  II,  and  ID  (late  time). 


AFWL  TRAPPED  ELECTRON  DATA  BASE 
BURST  :  STARFISH  SATELLITE  :  TELSTAR 
DAY  :  237  HOUR  :  1  MIN  :  0  ^  :  98 

CHANNEL  NUMBER  :  4  TIME  WINDOW  :  8.00 


Figure*  10  (a)  through  (g)  show*  raw  data  from  Explorer  XV  over  the  time  period  of  the 
Russian  2  and  3  bursts,  at  times  corresponding  to  those  shown  in  Figure  3,  for  comparison. 
It  is  apparent  that  fewer  points  are  plotted  from  the  AFWL  data  base  in  F  gures  10  (a) 
through  (g)  than  were  plotted  by  Roberts  in  Figure  3  for  the  1.9  Mc~  electrons,  even  though 
the  satellite  and  time  periods  of  coverage  were  the  same.  These  discrepancies  arc  not  resolved 
in  this  study,  although  the  previously  alluded-to  Grounding  procedure  may  account  for  the 
differences. 

The  temporal  progression  of  flux  after  one  burst  from  one  satellite  (Telstar),  at  a  typical 
L-value  of  1.9,  is  shown  in  Figures  11  (a)  through  fh).  Again,  typical  gaps  in  coverage  of 
Regions  I,  II,  and  111  are  readily  apparent  over  the  90  days  shown. 

As  previously  noted,  only  fundamental  and  two  higher  eignemodes  (/o,  J\ ,  /•>)  are  used 
to  fit  the  satellite  data  (Appendix  C).  The  third  eigenmodc  (/*)  is  not  used  becai  so  Shulr,  and 
Lanzcrotti  (Ref.  45:103)  suggest  that  the  modes  greater  than  two  rapidly  vanish.  Initially,  the 
steady-state  function  (fog)  was  used  along  with  the  fundamental,  but  this  solution  was  found 
to  be  “competing”  with  the  fit  for  the  fundamental,  in  the  sense  that  the  linear  constants 
tended  to  be  roughly  equal  in  magnitude  and  tended  to  alternate  in  sign.  Shulz  (Ref.  47:22 
23)  indicates  that  the  /«,  and  fo  functions  should  resemble  each  other  in  shape.  Also,  IMitr.er 
(Ref.37)  has  subtracted  the  background  from  all  flux  data,  and  the  background  should 
approximate  the  steady-state,  fag,  solution.  For  these  reasons,  the  steady-stale  solution  was 
not  used  in  the  data  fit. 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM  )/SEC 


INTEGRAL  ELECTRON  FLUX  (FLF.CTRONS/CM2)/SEC 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM  )/SEC 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM2)/SEC 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM2)/SEC 


BURST  :  RUSSIAN  2  SAT  :  EXPLORER15/EXPLORER15 
TIME  :  29  (dy)  19  (hr)  22.09  (min)  29.807(days) 
L  :  1.90  Beq  :  .045  BCUT  :  0.479  TW  (days)  :  4.0 
Bmin  (Gauss)  :  0.001  Bmax  (Gauss)  :  1.000 


0.5 

X  =  (1-Beq/B)1/2 


Figure  10(e).  Raw  Explorer  15  Data  at  Times  Corresponding  to  Figaro  3. 


INTEGRAL  ELECTRON  FLUX  (FLECTRONS/CM2)/SEC 


BURST  :  RUSSIAN  2  SAT  :  EXPLORER15/EXPLORER15 
TIME  :  47  (dy)  19  (hr)  22.06  (min)  47.807(days) 
L  :  1.90  Beq  :  .045  BCUT  :  0.479  TW  (days)  :  2.0 


0.5 

X  =  (1-Beq/B)1/2 


Figure  10(f).  Raw  Explorer  15  Data  at  Times  Corresponding  to  Figure  3. 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CKT)/SEC 


BURST  :  RUSSIAN  2  SAT  :  EXPLORER15/EXPLORER15 
TIME  :  63  (dy)  19  (hr)  22.06  (min)  63.807(days) 
L  :  1.90  Beq  :  .045  BCUT  :  0.479  TW  (days)  :  10.0 
Bmin  (Gauss)  :  0.001  Bmax  (Gauss)  :  1.000 


LEGEND 

a  =  .500  MEV 


io4  L- 


X  =  (1-Beq/B)'/‘ 

Figure  10(g).  tlaw  Explorer  15  Date  at  Timeo  Correaponding  to  Figure  0 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM  )/SEC 


BURST  :  STARFISH  SAT  :  TELSTAR 
TIME  :  2  (dy)  0  (hr)  0.00  (min)  2.000(days) 
L  :  1.90  Beq  :  .045  BCUT  :  0.479  TW  (days)  :  1.0 
Bmin  (Gauss)  :  0.001  Bmax  (Gauss)  :  1.000 


LEGEND 

a  =  .220  MEV 
o  =  .320  MEV 
a  =  .420  MEV 


X  =  (1-Beq/B)'^ 

Figure  11(b).  Temporal  Progression  for  Raw  Telstar  Data  from  1  to  00  Dap  Post-Burst. 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM2)/SEC 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CKr)/ SEC 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CMySEC 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM^/SEC 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM2)/SEC 


r# 

BURST  :  STARFISH  SAT  :  TELSTAR 
TIME  :  60  (dy)  0  (hr)  0.00  (min)  60.000(days) 
L  :  1.90  Beq  :  .045  BCUT  :  0.479  TW  (days)  :  5.0 


0  0.5  1 

X  =  (l-8eq/B)1/2 


Figure  11(g).  Temporal  Progression  for  Raw  Telstar  Data  from  1  to  90  Days  Post-Burst. 


INTEGRAL  ELECTRON  FLUX  (ELECTRONS/CM2)/SEC 


BURST  :  STARFISH  SAT  :  TELSTAR 
TIME  :  90  (dy)  0  (hr)  0.00  (min)  90.000(days) 
L  :  1.90  Beq  :  .045  BCUT  :  0.479  TW  (days)  :  5.0 
Bmin  (Gauss)  :  0.001  Bmax  (Gauss)  :  1.000 


LEGEND 

a  =  .220  MEV 
o  =  .320  MEV 
a  =  .420  MEV 


X  =  (1-Beq/B) ' 

Figure  11(h).  Temporal  Progression  for  Raw  Telstar  Data  from  1  to  90  Day s  Post-Rurst. 


70 


V.  RESULTS  AND  DISCUSSION 


Figure*  12  and  13  show  the  typical  temporal  progression  of  raw  data  at  two  representative 
L  values,  L  =  2.3  and  L  =  2.4.  Included  on  those  Ggures  are  the  flux  curves  calculated  from 
the  amplitude  coefficients  and  decay  times  fitted  to  that  data  by  Elcctrofit.  /'s  previously 
noted,  the  curves  shown  on  any  given  figure  may  vary  slightly  in  time  from  the  individual 
data  points,  since  the  curves  represent  only  the  one  liine-inslant  corresponding  to  the  start 
time  plus  one-half  the  time  window  of  the  plot.  Figures  12  and  13  represent  maximum 
freedom  of  the  fit,  since  each  r  is  free  to  vary  independently  of  the  other  r’s. 

Figure  14  shows  the  same  temporal  data  progression  with  different  fitted  curve:;.  The  curves 
in  this  figure  were  calculated  under  the  restriction  that  the  ratios  of  the  r’s  'verc  fixed  at 
the  theoretical  values  of  1:1/9:1/25;  hence,  only  T\  was  really  varying  in  a  random  sence. 

Figure  15  shows  similar  data  progression  and  fitted  curves  where  the  r’s  were  allowed  to 
freely  vary,  subject  to  the  constraint  that  ro  >  u  >  r->.  Also,  this  was  an  initial  run  with  500 
tries  for  convergence,  vice  96  tries.  In  addition,  the  computed  flux  was  used  in  the  iteration 
process  to  eliminate  non-physical  (negative)  solutions.  The  solution  was  restricted  to  being 
positive  at  four  different  times  and  at  five  different  equatorial  pitch  angles.  Alt  lough  this  is 
an  initial  run  with  only  a  small  amount  of  data  used,  it  indicates  that  the  recommendations 
in  Chapter  VI  should  be  pursued. 

Table  XVII  shows  the  fitted  linear  amplitude  coefficients  and  exponential  decay  times  (in 
days)  at  constant  Energy  and  Zcshell  for  a  representative  subset  of  the  data  examined  in 
this  study.  The  sheer  volume  of  resultant  fits  and  calculated  curves  made  it  impossible  to 
include  all  results  in  this  report.  The  values  shown  in  Tabic  XVII  were  computed  with  the 
r’s  allowed  to  vary  freely. 

Figures  12  and  13  are  considered  representative  of  the  results  of  the  fitting  model.  It  is 
apparent  from  examination  of  Figures  12  and  13  that  the  pitch-angle  diffusion  theory 
and  flux  computation  method  developed  in  this  report  are  consistent  with  experimental 
measurements.  Higher  eigenmodcs  are  readily  discernible  at  the  lowest  energy  in  the  Figures 
or  the  first  five  days,  and  are  not  visible  at  later  times.  These  higher  eigenmodcs  are  not 
discernible  in  the  plots  at  the  higher  energies;  Table  XVII  shows  that  these  modes  are 
present,  but  at  much  lower  levels,  in  some  cases.  The  calculated  curves  show  i  reasonable 
fit  to  the  data  even  at  late  times,  consistent  with  the  assumption  of  exponential  decay. 

Examination  of  Figures  14  and  15  shows  that  some  inconsistencies  remain  to  bo  resolved  in 
the  use  of  this  fitting  method.  Thu  caluclalcd  flux  curves  for  the  low-energy  data  points  also 
appear  to  be  slightly  lower  than  optimum,  while  the  high-energy  data  appear:  to  be  fitted 
better  by  the  model.  This  may  be  partially  explained  by  the  existence  of  many  lower  data 
points  in  the  vicinity  of  the  loss  cone  which  tend  to  “drag  down"  the  overall  fit.  Any  change 
in  the  cutoff  value  significantly  affects  the  fit  by  shifting  the  “vertical"  portion  of  the  curve 
toward  or  away  from  those  points.  However,  this  cannot  account  for  the  entire  difference 
of  the  curves  from  the  data.  The  energy  dependence  may  be  treated  less  Ilian  optimally 
in  the  choice  of  functional  solution  made  in  this  study.  If  the  chosen  function  il  solution  is 
not  the  correct  one,  il  may  be  nevertheless  close  enough  to  correct  to  give  a  r  rasonable  fit 
at  high  energy,  but  be  significantly  in  error  at  lower  energy.  A  different  functional  solution 
may  fit  the  lower-energy  groups  better  and  yet  still  lit  the  high-energy  elect  ons  as  well. 
Recommendation  for  the  use  of  other  functions  is  made  in  Chapter  VI. 
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Figure  12(a).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  f. 
Russian  3  Burst  ( r„  freely  varying). 
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Figure  12(e).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L 
Russian  3  Burst  (r„  freely  varying). 
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Figure  12(d).  Temporal  Progression  of  Raw  Data  r 
Russian  3  Barst  (r„  freely  varying). 
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Figaro  12(g).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Corn's  at  L 
Russian  .1  Hurst  (r„  freely  varying). 


BURST  :  RUSSIAN  3  SAT  :  EXPL0RER15  [045] 
TIME  :  25  (dy)  0  (hr)  0.00  (min)  25.000(days) 
L  :  2.30  Beq  :  .025  BCUT  :  0.506  TW  (days)  :  5.C 
Bmin  (Gauss)  :  0.001  Bmax  (Gauss)  :  1.000 


LEGEND 

a  =  .500  MEV 
°  =  .500  MEV 


X  =  (1-Beq/B) ' 


Figure  12(h).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L 
Russian  3  Burst  (r„  freely  varying). 
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Figure  13(d).  Temporal  Progression  of  Raw  Data  and  Fitted  Flnx  Carres  at  L 
Russian  3  Burst  (r„  freely  varying). 
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Figure  13(e).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curves  at  L  =>  2 A 
Russian  3  Burst  (r„  freely  varying). 
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Figure  13(f).  Temporal  Progression  of  Raw  Data  and  Fitted  Flax  Curves  at  It 
Russian  3  Burst  (r„  freely  varying). 
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Figure  14(b).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Carres  at  L 
Russian  3  Burst  (r„  fixed  1:1/9:1/25). 
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Figure  14(e).  Temporal  Progremion  of  Raw  Data  and  Fitted  Flux  Currw  at  L 
Roman  3  Bunt  (r»  fixed  1:1/9:1/25). 
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Fifare  14(f).  Temporal  PTogrewion  of  Raw  Data  and  Fitted  Flox  Carres  at  L 
Roman  3  Burst  (rn  fixed  1:1/0:1/25). 
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Figure  15(a).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Currrs  at  A 
Rumian  3  Burst  (rn  freely  Tarying  with  physical  contraints). 
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Fiffnre  15(b).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curres  at  L 
Russian  3  Oartt  (r„  freely  varying  with  physical  contraints). 
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Figv*  15(c).  Temporal  Progression  of  Raw  Data  and  Fitted  Flax  Curves  at  L 
Russian  3  Burst  (r„  freely  varying  with  physical  contrainUV 
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Figure  15(d).  Temporal  Progression  of  R3W  Data  and  Fitted  Flux  Carres  at  l, 
Russian  3  Durst  (r„  freely  varying  wit  h  physical  contraints). 
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Figure  15(h).  Temporal  Progression  of  Raw  Data  and  Fitted  Flax  Carres  at  L  =  2.4  for 
Russian  3  Burst  (r*  freely  varying  with  physical  contraints). 
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Figure  15(j).  Temporal  Progression  of  Row  Data  and  Fitted  Flax  C urm  at  L 
Russian  3  Burst  (rn  freely  varying  with  physical  co attaints). 
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Figure  15(m).  Temporal  Progression  of  Raw  Data  and  Fitted  Flux  Curre*  at  L  —  2.4  fo 
Russian  3  Burst  (r„  freely  raryinp  with  physical  contraints). 
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3.910  6  0033£>04  2  4441E«-03  1  910&t*C4  12  94  1  438  0  52  0  2716  4  11E-03 

7  3281E-06  3  0865E+06  3  5964EK55  3  1  0  2  923  2.28  0  2604  3  80E-03 
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(First  line  of  pair  is  initial  iteration;  second  line  is  result.)  (continued) 
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(First  line  of  p-.ur  is  initial  iteration;  second  line  is  result.) 


The  abort  time-spans  covered  by  most  of  the  data  sets  in  the  Trapped  Electron  Data  Base 
have  not  allowed  adequate  determination  of  P).  The  values  of  presented  in  Table  XVII 
are  suspect  for  this  reason.  However,  more  confidence  is  held  in  the  values  >f  r,  and  r>> 
computed  by  this  study.  In  most  data  sets  studied,  the  time  period  of  the  data  was  long 
enough  compared  to  the  values  of  rt  and  r>  to  compute  reasonable  values.  More  consistency 
between  data  sets  was  found  for  r(  and  r3  than  was  found  for  r0.  Of  course,  then  rn  may  be 
simply  computed  by  multiplying  r.  by  the  theoretical  ratio  corresponding  to  t tie  solution. 


112 


VL  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

The  pitch-angle  diffusion  theory  and  Sox  calculation  method  developed  in  this  study 
are  reasonably  consistent  with  experimental  data  in  most  cases.  The  method  shows  promise 
as  a  way  of  calculating  improved  decay  times  over  a  broad  region  of  the  magnetosphere, 
particularly  for  higher  energy  electrons.  Additionally,  a  comprehensive  summary  table  of 
the  best  available  literature  values  for  r0  has  been  prepared  (Table  XV).  A  complete  table 
covering  such  a  broad  range  of  energy  and  L  has  not  heretofore  been  compiled. 

Time  and  computer  limitations  have  precluded  examination  of  the  entire  data  base. 
However,  the  overall  efficiency  of  the  model  developed  in  this  study  has  been  demonstrated.  It 
is  expected  that  improved  predictions  of  fluxes  and  improved  inputs  to  the  AFW',  SPECTER 
Codes  (Ref.  7-9)  will  be  available  as  a  result  of  this  work.  This  will  enhance  l  ie  capability 
for  calculation  of  satellite  operational  environments,  and  will  improve  the  abili  .y  to  predict 
satellite  survivability  and  vulnerability. 


Recommendations 

The  following  recommendations  for  further  study  are  proposed  to  improve  the  model 
developed  in  this  study,  and  to  resolve  the  problem  of  non-physical  solutions  for  some  data 
sets.  At  the  direction  of  this  author,  several  of  these  recommendations  are  presently  being 
implemented  at  the  Air  Force  Weapons  Laboratory. 

(1)  Incorporate  an  estimation  of  the  error  in  the  magnetic  field  value,  D ,  of 
individual  data  points.  This  corresponds  to  a  “horizontal"  error  bar  on  the 
plots  such  as  Figure  12.  It  is  known  that  errors  exist  in  llie  compulation  of  the 
B  value  associated  with  each  data  point,  and  this  error  would  prove  significant 
in  the  region  of  the  loss  cone.  Improved  weighting  or  the  data  points  with  the 
combined  errors  should  improve  the  fit. 

(2)  Utilise  the  L-rounding  procedure  of  Chapter  IV  to  add  data  points  to  each  L- 
value  considered.  This  should  aid  in  filling  gaps  in  the  data,  whici  necessarily 
will  improve  the  fit. 

(3)  Discard  the  highest  eigenmode  if  the  fitted  value  of  its  decay  time  is  less  that 
0.25  day,  and  fit  the  data  with  only  the  fundamental  and  one  higher  mode. 
A  value  of  0.25  day  is  not  reasonably  fitted  within  the  errors  in  the  Trapped 
Electron  Data  Base.  This  mode  would  be  essentially  lost  within  thi  first  day  in 
any  case. 

(4)  Perform  an  “energy-scaling"  of  the  data  to  combine  data  from  different  satel¬ 

lites.  This  involves  assuming  some  form  of  energy  spectrum,  such  as 
including  an  estimate  of  the  error  in  the  assumed  value  of  F,q.  The  new  flux 
would  be  simply  the  old  flux  multiplied  by  a  factor  If  the 

error  in  the  old  flux  were  50  percent  (a  value  consistent  with  the  Data  Base) 
and  the  error  in  Eo  were  25  percent,  and  if  £o  «  1.1  (a  fission  spectrum),  then 
the  resultant  error  in  the  shifted  flux  would  be  less  than  1  percent  different 
from  that  of  the  old  flux,  when  shifting  from  080  KeV  to  500  KeV.  Thus,  it 
should  be  possible  to  combine  data  sets  and  fill  the  gaps  in  data  to  achieve  a 
better  fit.  However,  the  limitation  of  the  assumed  spectrum  still  exists. 
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(S)  Utilize  only  late-time  data  where  it  is  available,  and  (it  only  the  fundamental 
mode  to  it,  under  the  assumption  that  higher  modes  will  have  decayed  away. 
Then,  with  the  fundamental  mode  constrained,  go  back  to  early  time  data  and 
attempt  to  fit  higher  modes,  if  any.  This  recommendation  is  of  limited  value 
because  only  two  data  sets  cover  a  significant  period  of  time:  Slarfisb/Telstar 
and  Russian  3/Explorer  5.  However,  it  may  provide  a  better  value  of  rn  than  is 
presently  available  at  some  L  values. 

(8)  Revise  the  fitting  functions  to  other  Bessel  function  solutions  than  the  J-\fi 
solution  and  compare  results.  This  was  performed  only  on  one  oth'T  functional 
solution,  Jo,  with  only  one  small  data  set,  and  the  results  showed  no  improve¬ 
ment  over  the  present  study.  However,  to  be  certain  that  the  au  .hor  has  not 
simply  chosen  the  wrong  functional  form,  more  functional  forms  must  be  in¬ 
vestigated. 

(7)  In  the  fitting  program,  constrain  the  directional  flux  to  physical  reality  at  each 
iteration  (i.e.,  compute  flux  at  each  iteration  and  reject  non-physical  solutions 
at  each  iteration,  rather  than  after  converging  on  a  final  solution).  This  should 
converge  to  a  physically  real  solution  if  it  converges  at  all;  however,  significant 
increases  in  computer  processing  time  may  accrue. 

(8)  Increase  the  number  of  tries  for  a  solution  from  96  to  200  and  reduce  the  number 
of  outer  loops  from  32  to  10  or  IS.  This  may  improve  the  convergence  of  the 
least  squares  error,  if  non-physical  solutions  are  also  rejected  at  each  iteration. 

(9)  Compare  the  background  subtracted  from  the  flux  with  the  latest  available 
natural  environment,  to  ascertain  the  correctness  of  the  background  subtrac¬ 
tion.  If  significant  differences  are  found,  the  background  should  be  re-added  to 
the  flux  and  the  natural  environment  subtracted. 

(10)  A  suggestion  proposed  by  Professor  D.G.  Sbankland,  Air  Force  Institute  of 
Techuoiogy ,  involves  a  significant  modification  to  the  fitting  method.  This 
method  involves  an  interpolating  function,  or  measure  of  ‘roughness*  of  the 
fit,  which  may  be  made  arbitrarily  smooth.  The  advantage  of  this  method  is 
that  no  limiting  assumptions  need  be  made  about  tbe  data.  The  chief  disad¬ 
vantage  is  that  the  method  is  a  strictly  mathematical  method,  with  no  physical 
constraints.  Although  the  method  would  involve  extensive  reprogramming,  it 
promises  worthwhile  information  about  the  flux  curves,  and  should  be  tried. 
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APPENDIX  A 


Demonstration  of  Zero  Diffusion  Current 
at  the  Magnetic  Equator  for  Half-Integer 
Bessel  function  Solution 

To  demonstrate  sero  diffusion  current,  one  must  show  that 

limt_oZ)M  fn(z)  =  0 

for  the  (7  =  0  case.  We  Know  that  from  (31): 

J-\rM  =  >/*/(»*)  co*(2) 

and  from  (1ft): 

bMM  a*  bMMt 

and  from  (18): 


Hence,  the  deriratire  of  the  eigenfunction: 


M  -  -[^®EH][(S„ * „£]  + „ '  0 


Finally,  the  quantity  which  must  go  to  sero  is: 


The  term  in  square  brackets  is  a  finite  constant  for  any  n,  and  clearly  the  sine  goes  to  sero 
as  z  goes  to  sero.  Hence,  the  limit  is  sero. 


APPENDIX  B 


Derivation  of  Eigenfunction  Derivatives  7^  *nd  tfn 
To  find  the  derivative  of  the  steady  state  function,  one  starts  with  (Si): 

?oc= 

„_0  A» 

If  the  X*  are  assumed  (from  (19))  independent  of  x,  then  the  series  can  simply  be 
differentiated  term-by-term,  since  the  definite  integral  is  simply  a  constant  for  each  value 

7'-,*>-sLix!rs,-M^ 

”£*5 

n— 0  A"  JC 

To  find  the  derivative  of  the  eigenfunction  one  starts  with  (18): 


Using  the  product  rule  for  differentiation: 


(sm[(7»/2£e)(((l  -  v-)T0)fl  -  7C(1  -  vn/4))})] 

49~  j  9 


fo  —  -2* oo« 


-(*/»») 


W  [n/27^c(V^)[n/i  -  ^(«/^)][to  -  (To  -  r,)yV<jyj 
(s«n  |()r/(27e))(((»  -  V2)T0)/ 2  -  K (l  -  y'  /4))j)  dy 

/,  »  +2)ra,«-<,/T*)^[_^‘7ii[v/272e(3Jr/^)[>/l  -  y2(«/ftj] [tb  -  (To  -  7\)y*'4]y] 

(sin  [(3»/{2^))(((t  -  y-)r0)/2  -  Ar(l  -  y'  '4))1)  dy 

[N/27£(S3r/&)  \/l  -  y2(B/Bo)][r0  -  (T0  -  r,)yV«Jyj 
(sin  [(5»/(2Z,))(((l  ~  V-)To)/2  -  k[\  -  y'  /4))]j  dy 

h  =  *-2*a*e~i,'T')  '  \s/W^!2t)\\/l  -  y2(J?/Bo)]  [r0  -  (To  -  r,  )y»/4]y] 

(sin  |(7ir/(2Ze))(((l  -  y2)T0)/2  -  k(  1  -  ./  '“))])  dy 
Then  fitting  constants  are  really  ructions  of  energy  and  /y-shcll. 


/•»-> v/si/s 
h  =  -21^<l2*-<,/r,,  / 

*9—9* 


r,~y/lU/Bj 


The  ratios  of  ryj,  fj,  r2,  and  r;j  should  be  (as  computed  From  the  eigen  value  ratio  of  the 
zeroes  of  the  function): 


9f| 

»\i  ft*  25nj 
fb  w  49  rj 


ri  =  (l/0)n> 

or  alternatively  rj  =  (l/25)n» 
r3  =  (l/49)n, 


APPENDIX  D 


The  General  Solution  of  the  Integral-Energy 
Omnidirectional  Flux  Equation  for  Bessel  Function* 
of  Arbitrary  Order,  a*  a  Function  of  Sigma 

Substitution  of  (18),  (28),  and  (28)  into  (27)  completes  the  general  expre  ssion  for  the 
differential  energy  omnidirectional  flux,  J ,  with  Bessel  functions  of  order 

a  -  1 


where  a  is  some  number  less  than  two: 


•^(CirXir^-a-r") 
<T4«l»-ier . )] 

•vT'-yWM  r(vKsr2)) 

— M-(s) — ) 

♦(;)" . 4-criw-;X:n] 

•\A~-  »*(«/«,)  r(y)d(r))J 

The  analogy  between  the  above  equation  and  Equation  (43)  is  seen  if: 

jf  r'S/ih,.  dK 


and  if: 


[£" r*a„(/?)  '  {c  -  f/(r„))A  =  J~ p1  dK  A„,(K) 

=  f  p-  a„(/f)e-",%,,?,)  dE 
Jr 

Hence,  (AS)  is  simply  the  energy- integrated  form  of  (.17). 
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